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Abstract 

We review the status of the domain wall fermion approach to construct chiral 
gauge theories on the lattice. 

1 Introduction 

It is a natural question to ask, whether there exists a regularization of a chiral gauge theory 
beyond perturbation theory. This addresses the pure existence of the Standard Model of elec- 
troweak interactions. The astonishing answer to the above question is that a non-perturbative 
regularization has so far not been found. Even worse, nogo theorems Jl], Q seem to make it 
impossible to even find such a regularization. Of course, the tremendous success of the descrip- 
tion of the electroweak interactions in perturbation theory makes the Standard Model a very 
well tested theory in particle physics. However, this leaves us with the unsatisfactory situation 
of the Standard Model being good "for all practical purposes" . At this level it is certainly not 
a well founded theory as e.g. QCD. 

The conventional point of view today is to regard the Standard Model as only an effective 
theory to describe low energy physics. In this low energy limit perturbation theory works 
very well. However, due to triviality our theoretical description of the electroweak interactions 
is inherently incomplete. At some energy scale which depends on the Higgs boson mass the 
theory has to break down giving room for some -yet unknown- new physics. This new physics 
in some sense regulates the minimal Standard Model. Therefore, finding a non-perturbative 
regulator for a chiral gauge theory might give some hint about the theory "beyond the Standard 
Model" . A non-perturbative chiral regulator would also be important to understand, whether 
a construction of an asymptotically free chiral gauge theory is possible |3], [| . Not surprisingly 
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then, the attempts of a non-perturbative definition of a chiral gauge theory have been numerous 
and in particular the search for lattice chiral fermions has been an intensive area of research in 
the last few years. 

The outcome of all attempts to put the Standard Model on the lattice is, however, so far 
either negative or no convincing evidence of their success could be given. The stumbling block in 
constructing chiral gauge theories on the lattice are the famous nogo theorems by Karsten and 
Smit UU and by Nielsen and Ninomiya Q, see also ||. The assumptions of the nogo theorems 
are surprisingly mild. They require hermiticity, locality, and translational invariance of the 
underlying interaction. With these assumptions there will always be a clash between chiral 
gauge invariance and the possibility of a single Weyl fermion on the lattice. To each fermion 
with given chirality there will always be a doubler mode with the opposite handedness. This 
proliferation of fermions is deeply connected to the anomaly structure on the lattice. Either we 
will have chiral gauge invariance but then the doubler fermions appear and cancel the anomaly, 
or we lift the doublers by some mechanism and destroy chiral gauge invariance. It is worth 
to note that these theorems are claimed to hold also without refering to the lattice at all 0. 
This is reminiscent of Adler's demonstration that it is impossible to find a gauge invariant 
regularization such that the axial current remains conserved It emphasizes the fact, that 
the question of regularizing chiral fermions non-perturbatively is not a problem for the lattice 
alone but one of principle. 

To circumvent the nogo theorems one has then to violate some of the above assumptions at 
one stage. Although it is acceptable that at some intermediate step one of the above properties 
is missing, one clearly would like to recover them in the continuum theory as they are important 
ingredients for a sensible theory. There are numerous ways trying to find such violations and 
still recover the electroweak interaction in the continuum. For an overview I refer to several 
review talks in the annual Lattice proceedings j|, |7|, |[ || [UJ . In particular I would like to draw 
the reader's attention to the proceedings of the Rome workshop [O which contains virtually all 
approaches at that time. In [12| many of these attempts, the way how they tried to circumvent 
the nogo theorems and reasons why they did not work, are summarized. 

One attempt of finding a method to simulate chiral fermions on the lattice was not covered 
in the Rome workshop. This was Kaplan's idea [D| of generating a chiral mode on a domain 
wall that is induced by a mass defect along an extra dimension. Imagine a scenario where we 
actually live in a 5-dimensional world. Assume that along the fifth dimension there is a defect 
generated by a fermion mass term that only depends on the fifth direction and has the form of 
a soliton. This will generate a domain wall and on this 4-dimensional wall we will find a chiral 



zeromode traveling along it [[L4L 15 



In |Tj| Callen and Harvey studied general anomaly descent relations jl6] connecting fermion 
zeromodes on a domain wall and anomalies in even with parity anomalies in odd dimensions. 
In particular they demonstrated that the odd dimensional theory is anomaly free. Performing a 
Goldstone-Wilczek |17]] type of calculation they showed that the anomaly in the gauge current 
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due to the chiral zeromode on the domain wall is exactly cancelled by the divergence of the 
Chern-Simons current generated by the heavy fermions in the fifth dimension off the wall. That 
the Chern-Simons current is not divergence free is due to the fact that the mass defect changes 
sign when crossing the domain wall. 

One interesting aspect of the Callan-Harvey analysis is, that it can be taken over to the 
lattice. Indeed, it has been pointed out already in 113] that this might lead to a realization of 
chiral lattice fermions [} 

Independently of []TB| Kaplan suggested to use the domain wall setup to construct a chiral 
gauge theory on the lattice | 13fl . The hope is to produce a chiral theory on the 4-dimensional 
domain wall. Using the -in lattice QCD successful- Wilson mechanism, it can be expected that 
the dangerous doubler modes are decoupled. The difference to other approaches is that the 
5-dimensional theory one would start with is vectorlike from the very beginning thus avoiding 
problems with the nogo theorems. Through the Callan-Harvey mechanism the gauge anomaly 
would be reproduced on the domain wall while at the same time the 5-dimensional system 
remains anomaly free due to the mechanism of anomaly cancellation. Another attractive prop- 
erty of the domain wall model is, that also the fermion number current is anomalous on the 
4-dimensional wall [JT^, [T9fl , although, it is again conserved in the 5-dimensional system. 

This promising picture evokes immediately some questions. Staying in a situation where 
the fermions are free or are at most only coupled to external gauge fields, we have a single Weyl 
fermion on the wall. Then how did we circumvent the nogo theorems? A first answer is that 
we violated translational invariance. But this is only true in the 5-dimensional system. With 
respect to the 4-dimensional domain wall the system is translational invariant and still we have 
the single chiral fermion and the anomalies work out correctly. 

The answer to this seemingly paradox can be given, if one looks at the finite lattice situation. 
There we have to introduce some kind of boundary condition in the extra direction. This leads 
-unavoidingly- to a second domain wall on which an additional zeromode appears. This mode 
has the opposite chirality as the one on the original wall. Both modes have an exponentially 
small overlap fll3 , [201 21]. Therefore we are in full accordance with the nogo theorems. We 
indeed have the mode spectrum as demanded by the theorems. Due to the extra dimension 
we have just separated them such that they do not communicate. The whole setup is gauge 
invariant, but we created a loophole for the charges such that they can escape in the extra 
dimension leaving a chiral theory on the domain wall. For the infinite lattice the second 
domain wall is sent to infinity. That this picture is indeed correct, as long as we do not have 



dynamical gauge fields, could be demonstrated in 20, 22 



1 There a domain wall system was studied in the context of a stacking fault in a PbTe-type crystal. Using 
the argumentation of Callan and Harvey the authors of ref. Jl8[ discussed the apparent appearance of the parity 
anomaly. Using techniques similar to Callan and Harvey they showed that the complete 3-dimensional condensed 
matter system is anomaly free. They proposed that the domain wall model may be used for solving the chiral 
fermion problem in lattice gauge theory. In particular they suggest to introduce the Wilson mechanism for 
lifting the doublers. However, they then gave up the domain wall model idea. 
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There remains, of course, the crucial point what will happen when dynamical gauge fields 
are added. There have been two main lines of thought. The first one is to add 5-dimensional 
gauge fields giving the gauge fields in the fifth direction a coupling strength different from the 



4-dimensional fields [[L3| . A number of papers investigated this scenario |23|, [24], g^, [19], gg, |27 
The conclusion was that this setup would not lead to the desired chiral gauge theory. The main 
reason is that one would end up in the so-called layered phase where the system reduces to a 
standard 4-dimensional lattice system of Wilson fermions and no zeromodes appear. Hence the 
model becomes vectorlike and can not be used for a construction of a chiral lattice theory. 



The alternative approach pq , |25| , |29|| is to leave the gauge fields strictly 4-dimensional and 
to confine it to only one of the walls. In this way only one of the domain walls would be 
gauged and we are left with a region in the extra direction which is gauged and its ungauged 
complement. For reasons that will become clear later we call the gauged region the "waveguide" 
region. At the boundary of the waveguide to the ungauged part in the fifth direction gauge 
invariance is broken. This leads to a situation where one studies either a gauge non-invariant 
model or a model which is gauge invariant but contains scalar fields at the boundaries of the 
waveguide. The scalar (or Stiickelberg) fields serve to restore gauge invariance in much the same 
manner as the Standard Model is made gauge invariant. Adding a scalar field that couples to 
the fermions suggests, to also introduce a Yukawa-coupling y. It can therefore be expected 
that fermion masses are generated via the Higgs mechanism through the Yukawa-coupling of 
the fermions to the scalar vacuum expectation value. These fermion modes will live at the 
waveguide boundary. The crucial question is, whether these states can be decoupled such that 
we are left with only the gauged and ungauged domain wall zeromodes separated from each 
other in the fifth dimension. 

One may argue that such a scenario is condemned to fail from the very beginning. The reason 
being that approaching the phase transition from a spontaneously broken to a symmetric phase 
the vacuum expectation value v approaches zero giving rise to small fermion masses which may 
be interpreted as mirror fermions. In particular these masses should be zero in the symmetric 
phase. Thus they would not decouple. However, investigations of Yukawa models on the lattice 
revealed a surprisingly rich phase structure at large values of the Yukawa-coupling [RD . In 
particular a so-called strong coupling region could be identified |H[ [32], p3| . In this region the 
Yukawa-coupling becomes so strong that the fermions combine with the scalar fields to form 
massive bound states with masses at the cut-off. If such a region would exist also in the domain 
wall model, the fermions at the waveguide boundary could be made heavy and would decouple. 
We then would be left with a single chiral zeromode bound to the domain wall. 

Another argumentation which would make the domain wall model fail is that the charge 
due to this chiral zeromode will flow off the wall transported by the Chern-Simons current and 
would have to be deposited somewhere. The most easy scenario to imagine is that there are 
additional fermion zeromodes that induce a current absorbing this charge pi ]. However, due 
to our waveguide setup, the gauge fields have to stop somewhere in the fifth dimension. At this 



4 



point we will have a change of the gauge field that can induce a Wess-Zumino current. This 
current can give rise of cancelling the anomaly. 

From these considerations it is clear, that the waveguide model can not be ruled out by 
simple arguments. The question whether the domain wall approach would be successful or not 
becomes a dynamical one and a numerical simulation becomes necessary. Such a simulation for 



a 3-dimensional setup was performed in There the weak Yukawa- coupling region, where 
the fermion masses are proportional to v and therefore vanish in the symmetric phase, could 
be established. For large values of the Yukawa-coupling the results of this simulation were 
compared with Yukawa models that have been investigated earlier and that do or do not show 
a strong coupling behaviour. Strong evidence was collected that in the domain wall model 
a phase with the strong coupling behaviour, as described above, does not exist. This was 
confirmed by the results of a strong coupling expansion which revealed a new phase with weak 
coupling behaviour at large y. It has to be concluded therefore that in the domain wall model 
physics is governed by a weak coupling behaviour at all values of y thus not leading to a chiral 
theory. 

In this review we will in the following discuss the domain wall model in more detail. 
Throughout this paper the fermions are thought to be taken in an anomaly free represen- 
tation. After introducing the fermion doubling problem and how it may be solved by the 
domain wall model we will continue to give an explicit demonstration of the Callan-Harvey 
picture of charge flow on the lattice. We then present the model with the scalar fields added. 
We discuss a large y expansion in this model and the results of the numerical simulations. We 
will present the suggestion that the extra dimension can be held strictly infinite. We end with 
a possible prospect of the domain wall model for simulating QCD and give a summary and 
conclusions. 



2 Free fermions, doubling and a possible escape 

In this section the problem of fermion doubling on the lattice |35| is illustrated and the idea 



of domain wall fermions introduced. The discussion will be for two dimensions only and the 
Hamiltonian formalism will be used. The results obtained in this simple setting show already 
the essential features and are readily generalized to arbitrary dimensions. In this section we will 
indicate the lattice spacing a explicitly in the formulae. Otherwise a is set to one throughout 
the paper. We start our discussion with the 1-dimensional Hamilton operator for free fermions 
which -in the continuum- may be written as H c = — ffi(^ + mo). The operator (j) = <72<9/cte acts 
only on the x coordinate and the <r's are the usual Pauli matrices. H c describes a single fermion 
of mass m obeying the relativistic dispersion relation E 2 = p 2 + m^. A "naive" discretization 
of H c on the infinite lattice would be 

H = -o x [a 2 d x + m ] , (1) 
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where d x now denotes the lattice derivative 

9x = [5z,x+^ $x,x— fj\ 1 (2) 

a denotes the lattice spacing and \i is a displacement by a in a given direction on the lattice. 
In our 1-dimensional example \i — a. 

By Fourier transformation, the particle energies are easily obtained as 



E = ±^ sin 2 (ak) + m 2 , (3) 

with k lying in the Brillouin zone —it < ak < it. For ka <C n, the naive continuum limit, a — > 0, 
reproduces the above relativistic dispersion relation E 2 = k 2 + m 2 ,. However, for momenta 
ak ~ it, we find again E 2 = k 2 + TOq and obtain another fermion with mass mo- Moreover, 
while the fermion at the origin of the Brillouin zone at ak ~ is a right moving particle as 
can be seen from the group velocity dE/dk, the one appearing at the corner of the Brillouin 
zone at k « rr/a is a leftmover since the sin function changes its sign there. For massless 
fermions this amounts to find two chiral fermions with opposite chirality quite in contrast to 
the single Weyl fermion we would have started with in the continuum. This is the famous 
doubling phenomenon for lattice fermions. The fact that one finds opposite chirality fermions 
at the different corners of the Brillouin zone is a consequence of the "doubler symmetry" on 



the lattice |36]]. This symmetry transformation is represented by a matrix 

M = M 1 M 2 ; Mj = ia 3 aj . (4) 

If ipk is an eigenfunction of the Hamiltonian in momentum space, this transformation exchanges 
the corners of the Brillouin zone, ipk = Mip k+7T = —<r 3 ip k+7T . Thus, the mode at the opposite 
corner of the Brillouin zone has a flipped chirality. The doubler symmetry gives in general 
dimensions an equal number of left and right handed particles. As the representation of the 
corresponding doubler symmetry group is irreducible, doublers have to appear in the fermionic 
lattice spectrum as long as the doubler symmetry is unbroken. 

The solution of the doubling problem as proposed by Wilson H] is to break the doubler 
symmetry by adding a second derivative term to the Hamiltonian 



with 



The spectrum becomes 



E 
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H = -o\ [<J 2 d x + m - rAj , (5) 



A x — — [6 X)X +n + 5 XjX —n — 25 X)X ] . (6) 



1 / r \ 2 

— sin 2 (afc) + ( m H — (1 — cos(afc)) J . (7) 
a 2 \ a J 
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The mass of the fermions are now given in the a — > limit as 

T 

m = mo + 2— rijr (8) 
a 

with n-x = for the origin of the Brillouin zone at ak « and n n = 1 for the corner of the 
Brillouin zone at ak w 7T. Therefore, we obtain for the fermion at the origin of the Brillouin 
zone again the continuum dispersion relation but the fermions at the corners get an additional 
piece proportional to 1/a. They will become infinitely heavy as a — > and decouple from the 
spectrum. In general dimensions counts the number of it's at the different corners of the 
Brillouin zone. 

Unfortunately, there is a price to pay for getting rid of the extra unwanted fermions. Imagine, 
we add dynamical gauge fields to the Hamiltonian (|J) by making the lattice derivatives (@),(@) 
gauge covariant. The Wilson term clearly acts as a mass term so that even if we take the 
fermion mass mo to zero, we still loose chiral gauge invariance. Although chiral symmetry can 
be restored in QCD by carefully tuning the bare parameters to some critical value, this appears 
to be a disaster for the electroweak interactions. 

Indeed, for vectorlike gauge theories as QCD, weak coupling perturbation theory revealed 
the restoration of chiral symmetry for the Green functions in the continuum limit [0, |37j . This 
opened the road for lattice simulations of QCD because -though technically very challenging- no 
problem of principle remains. For a chiral gauge theory using Wilson's approach [^, [7], |38|, weak 
coupling perturbation theory |39| showed no success in obtaining the desired target theory and 



the doublers could not be decoupled. This pessimistic picture was strengthened by results from 
lattice simulations which showed that for weak coupling the doublers remain in the spectrum 



11, 40 1 . Also the hope of decoupling as a non-perturbative effect failed [41, |42j|. It seems that 



in order to obtain a chiral gauge theory one has to be more ingenious than in QCD. 
2.1 Domain wall fermions 

A possible escape from the doubling problem is due to D. Kaplan who suggests to send the 
dangerous modes into an extra dimension. It is known since a long time that there exists a 



chiral zeromode solution of the continuum Dirac equation in presence of a soliton |14[ . The idea 
to use such a kind of solution to construct chiral lattice fermions is to regard our 4-dimensional 
world as a domain wall embedded in 5 dimensions. The interface is induced through a mass 
defect in form of a soliton. To stay in our 2-dimensional setup of the previous section we will 
study a 3-dimensional system. Along the third extra dimension denoted by s, a mass defect is 
introduced through 

f \ J ~ m ° ; S ~* ~°° (n\ 
m(s) = < . (9) 

[ +m ; s — > +oo 

The concrete form -aside that it has to be monotonic- of the mass is not very important. It is 
easy to show that in this situation there exist solutions ^± which are energy eigenstates of the 
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continuum Hamiltonian 

H c = -a 1 [a 2 d(x) + a 3 d(s) + m(s)} , (10) 
V± = e ip * x <S>±(s)u± (11) 

with 

$± (s) = exp (± jf m{s')dA (12) 

and u± a chiral eigenstate, a 3 u± = ±u±. Only the function $-(s) corresponds to a normalizable 
solution. This solution describes a chiral zeromode traveling along the interface. It is bound 
to the wall and falls off exponentially with growing distance from the wall. 

Translating this model to the -infinite- lattice we choose a step function for the mass 

!— 1; s < —a 
0; s = . (13) 
+1; s > a 

The continuum derivatives are again replaced by finite lattice differences and the Hamiltonian 
reads now 

H = -o x [a 2 d x + a 3 d s + m(s)] . (14) 
Imposing an ansatz similar to the continuum solution 

^ ± = e lkx §±{s)u ± (15) 

we find again a normalizable solution 

$_( s ) = e"^ 1 , (16) 

with sinh /io = mo. This is a solution that is bound to the domain wall, it is chiral and describes 
plain waves along the direction of the domain wall. Unfortunately, there exists another solution 
which only appears on the lattice 

^ + = e ikx <5> + u + (17) 

with 

s + 00 = (-i)'M*) ■ (is) 

This solution describes again a massless fermion traveling along the domain wall, it is also 
bound to it and has opposite chirality as compared to the solution $_(s), eq . flllf ) . In other 
words it is a doubler fermion in the s-direction! Even worse, also the doubler mode in the 
x-direction from the corner of the Brillouin zone is a solution. Therefore the spectrum on the 
domain wall consists of two chiral fermions with positive chirality and their two doublers -one 
in the s and one in the x direction- with opposite chirality. 

What we have achieved so far is instead of removing the doubler modes we have created one 
in addition. But what if we now try to apply Wilson's trick adding a higher derivative term? 
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The 3-dimensional model we started with is vectorlike as in QCD and we saw that there the 
Wilson term is harmless. The Hamiltonian becomes 



H = -<r x [a 2 d x + a 3 d s + m(s) - r{A x + A,)] . (19) 

Let us set the Wilson coupling r = 1 for the moment and again try an exponential ansatz for 
the transverse wavefunction $ keeping plane wave solutions in the x-direction 

$±(s±a) = -k eff (s)$ ± (s) (20) 

with an effective momentum 

k eff {s) = m B{s) - 1 - F(k); F(k) = 1 - cos(a£;) . (21) 

In order to obtain normalizable solutions, the functions $ have to decrease with growing \s\ on 
both sides of the wall. This implies the following conditions to be fulfilled 

$+ : \k eff \ > 1 for s < 0, \k eff \ < 1 s > 
$_ : \k e ff\ < 1 for s < 0, \k e ff\ > 1 s > . 

One recognizes that only one of the solutions is normalizable namely $_ and that we have to 
throw away the $ + solution. Thus we are left with only one zeromode solution and we got 
rid of the doubler mode in the s-direction ([18]). Since the doubler mode in the x-direction is 
decoupled by the usual Wilson mechanism we remain with a single chiral fermion solution as in 
the continuum. This chiral zeromode is bound to the domain wall and exponentially decreases 
going away from the wall. The normalizability condition (^) requires 

< m - F(k) < 2 . (23) 

This means that for a given value of the domain wall height mo the chiral zeromode does 
not exist for all momenta but only up to some critical momentum k c . At this value the 
fermion ceases to be chiral and vanishes in a band of heavy modes. This behaviour is the 
important ingredient of the domain wall model. The chiral zeromode only exists as a low energy 
phenomenon. For large values of the momentum it becomes heavy and is not distinguishable 
from the doubler modes. It therefore does not appear again as a low momentum state at some 
other corner of the Brillouin zone. 

In PI| the zeromode solutions for general m and r values were computed. There the 3-di- 
mensional Dirac-operator was studied. Therefore the momentum space becomes 2-dimensional. 
The normalizability condition led to a very peculiar zeromode spectrum (see also fig.l in [pT]| ) . 
For m/r < no chiral fermions exist. Increasing m/r a chiral fermion is found with the critical 
momentum given by m = rF(k). Increasing m/r further the value of the critical momentum 
keeps growing, too. This continues until at m/r = 2 the maximal critical momentum is reached. 
For m/r > 2 the chiral zeromode at the origin of the Brillouin zone k = (0,0) is lost. Instead 

9 



new zeromodes appear at the corners of the Brillouin zone k = (tt, 0) and k = (0, 7r). They 
have the opposite chirality as the one at k rs 0. The critical momenta of these new zeromodes 
are in the interval determined by m = rF(k) for the lower and m = r(F(k) + 2) for the 
upper momentum interval boundary. Increasing m/r even further we loose these modes again 
for m/r > 4 but get a new zeromode at k — (tt, tt) with flipped chirality. Finally, for m/r > 6 
the chiral zeromodes disappear completely. Depending on the bare parameters of the domain 
wall model one can therefore generate a chiral zeromode spectrum as needed. 

2.2 The domain wall model on the finite lattice 

The interesting question is, of course, whether the above scenario can be realized on a finite 
lattice. There one has the additional complication that some kind of boundary condition has to 
be imposed which necessarily induces a second domain wall. In the following we let the lattice 
have finite extensions L s in the s and L in the x-direction and have s G [1, L s ], x G [1, L\. The 
domain wall mass will be chosen 

f -1 2 < s < \ 

mis) = sinhOo)0(s); 9(s) = 1+1 I f + 2<s<L s . (24) 

{ 8 = 1,^ + 1 

The finite lattice Hamiltonian including the Wilson term reads 

H = -a l [a 2 d x + a 3 d s + m{s)-r{A x + A s )} . (25) 

Due to the finite lattice extension, we have to specify some sort of boundary condition which 
will be chosen to be periodic in the s and antiperiodic in the x-direction. The Hamiltonian (|25|) 
can be reduced to only depend on s by imposing plane wave solutions in the x-direction. One 
obtains, setting a = 1, 

H = -a x [a 2 sm(k) + a 3 d s + m(s) +r(cos(k) - 1) - rA s )} , (26) 
where the momenta k are now discretized, 

k = l( n + \)> n = 0,...,L-l. (27) 

This choice of the momenta k corresponds to antiperiodic boundary conditions which will also 
be chosen in numerical investigations presented later. They become necessary in numerical 
simulations to avoid exact zeromodes which render the simulation algorithms impractical. 

In this form the energy eigenvalues and the corresponding eigenfunctions for a given value 
of k can be obtained by diagonalizing the Hamiltonian which is a 2L S (g> 2L S matrix numerically. 
In this way one gets the momentum dependent eigenvalues X±k and their corresponding wave- 
functions. The eigenvalues come in ± pairs. In fig. la the wavefunctions belonging to the lowest 
positive eigenvalue Ao and its negative partner are plotted. One clearly sees that the second 
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Figure 1: The zeromode wavefunctions as obtained from the finite lattice Hamiltonian fl26). 
The domain wall mass = 1, the Wilson coupling r = 1 and the lattice size in the x-direction 
is L = 32. The horizontal lines indicate the position of the domain walls. Fig. la is the domain 
wall model with an extension of L s = 100 in the extra dimension. Fig. lb is the variant of this 
model with open boundary conditions and an extension of L s = 50 only. Both figures show 
that the zeromode wavefunctions are sharply localized on the walls and fall off very rapidly 
with growing distance from the wall. 



domain wall generates an additional solution which is absent on the infinite lattice. The nice 
thing is that both solutions are located on each of their domain walls at s = 1 and s = L s /2 + l. 
They have opposite chirality and fall off exponentially going away from the domain wall having 
only an exponentially small overlap oc e _W)is . 

In fig. 2a the lowest and the next to lowest positive eigenvalues are plotted as a function of 
the lattice momenta k corresponding to antiperiodic boundary conditions. One observes the 
existence of the critical momentum k c . The lowest eigenvalue Ao follows the dispersion relation 
of a free lattice fermion sin(/c) up to the critical momentum where it vanishes in a band of larger 
eigenvalues. In fig. 2a also the next lowest positive eigenvalue Ai is plotted. For low momenta 
where Xo{k) = sin(fc) it stays at the order of the lattice cut-off exhibiting a large gap to the 
lowest eigenvalue. It stays almost constant as a function of the momenta until at k c 1 it 
combines with X . Fig. 2b shows that the mode corresponding to the lowest positive eigenvalue 
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really is a chiral fermion. The ratio build from the wavefunction ip belonging to Aq 



R = ZM* > (28) 

indicates the chirality of a mode. It is zero when ip is a chiral eigenstate and becomes non-zero 
if chirality is lost. The figure shows that this happens exactly at the place where the gap Ai — Ao 
in fig. 2a shrinks to zero which is at k c . 



3 The anomaly 

The fermion spectrum we have obtained in the previous section looks very promising. We have 
two Weyl fermions on the finite lattice living on each of the domain walls. They are separated in 
the extra dimension and have an only exponentially small overlap. Besides the chiral fermions 
we have the heavy doubler modes on the walls and additional heavy fermions off the wall with 
masses at the order of the cut-off. For free domain wall fermions the low energy physics on 
both walls appears therefore to describe a chiral theory. 

Let us assume that adding gauge fields will not destroy this picture. This is, of course, 
a strong assumption but, as we will see, it is at least justified for weak external gauge fields. 
Introducing dynamical gauge fields is much more problematic and will be discussed later. A 
(2-dimensional) physicist living on one of the walls knowing nothing about he extra dimension 
would eventually discover "his" Standard Model - the 2-dimensional chiral Schwinger model. 
He would also discover the chiral anomaly and wonder how he could regularize his chiral theory. 
If the domain wall model is able to succeed in producing a chiral gauge theory on a wall as 
the low energy limit of the 3-dimensional model, it should as a first step generate the correct 
anomaly structure for external gauge fields. Therefore the computation of the anomaly in the 
2-dimensional Schwinger model is a first crucial step in testing the domain wall approach. Since 
at the present stage the gauge fields are only external, we might select our physical wall by 
hand. It is therefore irrelevant that the gauge fields connect both chiral zeromodes. What we 
are aiming at, is to get a picture of the charge flow in the finite lattice model containing both 
walls. For the following calculations we use a variant of the domain wall model as introduced 
by Shamir H3| and follow H] to compute the axial charge in the Hamiltonian formalism from 



the spectrum flow. 

Shamir's variant is to not let the mass change its sign across the second domain wall but 
to cut it off completely and introduce open boundary conditions instead. The system can 
be regarded as being in a box with infinitely high walls at the ends of the extra dimension. 
The Hamiltonian analysis of this situation is very similar to the above discussion fl43| . In this 
setup the chiral zeromodes appear as surface modes on the walls with opposite chirality on the 
two borders. We show in fig. lb the zeromode spectrum of the Hamilton operator with open 
boundary conditions in the s-direction. As before we keep antiperiodic boundary conditions 
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Figure 2: The lowest and next to lowest positive eigenvalues are plotted in fig. 2a against the 
lattice momenta k = 7i(n + 0.5)/ L with n = 0, ...,31 corresponding to a lattice of 32 points 
in the x-direction with anti periodic boundary conditions. The parameters are the same as in 
fig. 1. The dotted line is the energy of a free massless fermion on the lattice, sin(fc). One nicely 
observes that the lowest eigenvalue Ao (full squares) follows the sin(fc) behaviour whereas the 
next mode Ai (open squares) is at the order of the cutoff and stays almost constant as a function 
of k. At a critical value of the momentum k c ~ 1 the gap Ai — Ao shrinks to zero. At this point 
the lowest mode departures from the sin behaviour. Fig. 2b shows that the mode corresponding 
to Ao is really chiral. Plotted is the ratio R =< ipip > / < ipuiip > for the eigenfunction of Ao- 
R is zero if ip is a chiral eigenfunction and deviates from zero if the chirality is lost. R becomes 
non-zero at exactly the value of k where the two modes in fig.2a combine, i.e. at k c . 
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in the x-direction and choose the same parameters as in fig. la. The figure shows the surface 
modes which are again separated in the extra dimension. They have only exponentially small 
overlaps and are both eigenstates of 03, the chiral operator. It is remarkable how similar they 
are in shape compared to the domain wall zeromodes depicted in fig. la. The boundary version 
of the domain wall model has the advantage that the lattice in the s-direction is only half that 
of the domain wall model. This might lead to an improvement in practice if one thinks of 
numerical simulations. This would be the case in particular when the domain wall fermions are 
going to be used for simulations of QCD, a prospect we will discuss later. 

To see that there really flows an anomalous current along the domain wall we apply a time 
dependent external field in the x-direction. We introduce the gauge covariant lattice derivatives 

d x (U) = l - [U X J X>X+ , - U:_^5 X ^ X ] (29) 

and 

^x{U) = - [U X;lx 6 XjX+lx + U*_^5 X ^ >X - 25x, x . (30) 

Here U denotes the usual compact representation of the gauge field. In this review we will for 
all discussions only use U(l) gauge fields. Then U is a phase, U Xjll = e^ iClx '^ and is related to 
the continuum gauge potential A^(x) by e tqA ^ x \ The gauge fields will be chosen to be only 
2-dimensional and identical on every s-slice. Although this would lead to an interaction of 
the surface modes when the gauge fields were fully dynamical, this is not a problem for only 
external gauge fields because in this case a domain wall can be singled out by hand. 
The Hamilton operator becomes 

H = -a x [a 2 d x (U) + a 3 d s + m(s) - r(A x (U) + A s )} . (31) 

If we think of H being written in terms of creation and annihilation operators t/m an d V>, it is 
gauge invariant under a gauge transformation 

U X)li -> g x U x ^g* x+fl ; ip x -> g x ip x , g G U(l) . (32) 

We will choose an external field that is constant in space and varies only in time. For adiabatic 
external fields we expect then from the continuum 



Q 3 (t) = fdt 1 f dxe^F^(x,t') 

Z7T JO J 

2-7T JO OV 



= -2q±-A x {t). (33) 
Here we identify the gauge potential A\(t) = qa(t), with q — 1 the gauge charge and Q3 denotes 



the (axial) charge. Due to the gauge invariance (|32D the Hamiltonian (31) depends only on 



the product of all the gauge links. Thus the gauge potential U = e ta becomes a constant 



14 



global phase factor acquired by a fermion traveling around the finite lattice system. The lattice 
Hamiltonian reads in momentum space 



H = —<7i [a 2 sin(fc — a) + a 3 d s + m(s) + r(cos(A; — a) — 1) — rA s )] . (34) 
This is the free Hamiltonian with momenta shifted by the phase a. 



A definition of the axial charge is quite arbitrary on the lattice ||45|| . Basically one has to 
fulfill that opposite charge has to be assigned to the two chiral modes. Defining the vacuum as 
the state with all negative energy levels filled, we can define a measure of the charge by 

Qs = y~~~~\ £( L - " 1 - 2s XA* ( 35 ) 

— 1 q,s 

where the sum goes over the one particle wavefunctions ip s ^ q above the vacuum. The definition 
of Q 3 in (|3"5D has the property to be one, if the mode is located exactly to the left domain wall 
and minus one if it is located exactly to the opposite wall. For heavy modes that are smeared 
over the lattice the contribution to Q3 is zero. In [44] the charge Q3 was measured. The result 



is shown in fig. 3 where we plot the charge Q3 as a function of the flux a. If we measure a in 
units of 2n/L the anomaly equation ( |3~3"D becomes Q 3 = —2a. This behaviour can indeed be 
observed in fig. 3. 

If the external field would really be adiabatic, the turning on of the field would be so slow 
that the surface modes on the walls can tunnel through the extra dimension, exchanging there 
roles. This would happen at a = 1/2. At this point we will have a change of the zeromode 
spectrum (see section 2.1). A filled level on one of the boundaries is just about to become the 
positive energy particle and an empty level is about to drop in the sea. We would observe a 
jump in the charge by —2 and the system would relax to the ground state at a = 1 resulting 
in a zero net charge. This is indicated in fig. 3 as the solid line. 

We saw, however, that the wavefunctions of the surface fermions are sharply peaked at the 
boundary, see fig. lb. Therefore their overlap is exponentially small resulting in a very tiny 
tunnel energy 6. For the situation in fig. 3, 5 ~ 10 -9 . Thus 1/5 is the largest timescale in the 
problem. Instead of tunneling of both surface states, a particle-hole pair is created and the 
charge keeps decreasing towards Q 3 = — 2 at a = 1. This is represented as the stars in fig. 3. 



3.1 The resolution of a puzzle 

For the reader who knows of the extra dimension there appears to be a puzzle from the results 
of the previous section. In every step of the above discussion the odd dimensional model was 
gauge invariant and vectorlike. We are also in full accordance with the nogo theorems in that 
there are an equal number of left and right handed fermions. The only trick to circumvent 
them is that these modes are separated in the extra dimension. In fact, the whole idea was 
to study a vector theory in order to be able to use the Wilson mechanism in close analogy to 
lattice QCD to decouple the doublers. 
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Figure 3: The charge Q3 as function of the "flux" a. The slope is ~ —2 as expected from the 
continuum anomaly equation. 

So, how can there be then anomalous charge creation? The resolution of this puzzle are the 
anomaly descent relations [[TjJ [16[ . Callan and Harvey pointed out that in an odd dimensional 



theory with a mass defect in form of a domain wall the heavy modes residing in the extra 
dimension off the domain wall do not entirely decouple but produce a Chern-Simons term 
in the effective action when they are integrated out. This action generates a Chern-Simons 
current which is very similar to a Hall-current in solid state Physics. It flows along the extra 
dimension and perpendicular to the gauge current induced by the external field. The Chern- 
Simons current is responsible for a charge transport in the extra dimension from one of the 
walls to the other. The strength of the Chern-Simons current is different on the two sides 
of the wall. Therefore the Chern-Simons current is not divergenceless. Indeed, as could be 
shown by Callan and Harvey, the divergence of this current is exactly equal in strength to the 
anomaly on the wall. This mechanism makes the even dimensional observer think that charge 
just "disappears" or is "generated" out of nowhere. The total 3-dimensional system is anomaly 
free and contains an equal number of left and right handed fermions just as demanded by the 
nogo theorems. 

With these remarks we leave now the Hamiltonian formalism to discuss the domain wall 
model and go to the euclidean action formulation in which simulations have to be performed in 
the end. We start with an explicit lattice computation of the Chern-Simons current, confirming 
the continuum picture of the anomaly cancellation by Callan and Harvey being realized on the 
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lattice, too. 



4 Euclidean lattice action for domain wall fermions 

The Hamilton formalism for free domain wall fermions discussed in the previous two sections 
led to some remarkable properties of the domain wall model in its wall-anti-wall realization on 
a finite lattice. There is a single Weyl fermion bound to the domain wall with its chiral partner 
living on the anti domain wall separated through the extra dimension. All other doubler modes 
are heavy and decoupled by means of the Wilson mechanism. From the spectral flow we saw that 
the charge on the domain wall very closely followed the continuum anomaly equation. These 
properties are certainly very essential building blocks should the model succeed in describing 
a chiral gauge theory on the lattice. Would it have failed at this point, there would have been 
no reason to continue with its investigation. The considerations of the last sections can be 
regarded as some necessary first checks every lattice chiral fermion proposal should pass in 
order to have a chance of being successful. 

Our aim is a non-perturbative understanding of a chiral lattice model. This will most 
probably involve a numerical simulation at one stage. Such simulations are done by evaluating 
the -euclidean- path integral 



with Sq the gauge field action and Sp the fermionic part. We will not give the explicit form of 
both actions. Various alternatives will be discussed in the next section when possible additions 
of dynamical gauge fields to the domain wall model are discussed. Since we are aiming at a 
euclidean lattice description, we switch now the language going from the Hamiltonian formalism 
to the path integral formulation of the model. We start by performing an explicit calculation 
of the Chern-Simons current on the lattice. We will in the following stay in the d — 2 + 1 
dimensional setting and will always have abelian gauge fields. 

In their continuum calculation, Callan and Harvey pointed out |15] that in a situation with 
a mass defect in odd dimensions there is a Chern-Simons current which is generated by the 
heavy fermions living far off the wall. This current has opposite sign on the two sides of the 
wall and its divergence is hence non-zero. In fact, Callan and Harvey computed the divergence 
of the Chern-Simons current and found that it exactly reproduces the coefficient of the anomaly 
generated by the chiral zeromode on the domain wall. It is this mechanism which gives hope 
that the domain wall model is able to provide the anomalies which are needed for a reproduction 
of the anomaly structure in the Standard Model. 

The Callan-Harvey analysis can be considered to be incomplete in two aspects, (z) The 
computation of Callan and Harvey has been done far off the wall and a computation also close 
to the wall would be desirable, (ii) The Chern-Simons coefficient is regulator dependent and it 
should be computed in different regularizations. In fact, the continuum calculation of Callan 
and Harvey was done without imposing a regularization at all. We will improve this situation 




(36) 
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by presenting an analysis of the Chern-Simons current that is valid also arbitrary close to the 
wall |56j and give an explicit lattice computation of the Chern-Simons current [O, |22|, ^5 



(i) The effective Chern-Simons action is obtained when the heavy fermion modes are in- 
tegrated out. The gauge variation of this action then leads to the Chern-Simons current the 
coefficient of which can be calculated in the low energy limit. This program has been gone 
through in the continuum by Callan and Harvey. However, their computation was done far off 
the domain wall where the mass defect can be considered as constant. It would be obviously 
desirable to have an analogous calculation also arbitrary close to the wall Q. 



Such a continuum calculation has recently been given in [|||. Choosing the mass defect in 
the form 

m(s) = m tanh(m s) (37) 

we search for solutions of the eigenvalue equation in Minkowski space <j\ {ia^d/dx^ + m(s)} ip\ = 
Xipx, fi = 1,2,3. The problem of solving this equation for the 3-dimensional free domain wall 
system can be transformed to that of a quantum mechanical scattering problem with a modified 
Poschl- Teller potential which can be solved exactly. Thus we not only get the chiral zeromode 
solution discussed in section 2, but also all the heavy fermion solutions. From the wavefunctions 
ipx and the eigenvalues A the propagator G(z, z'), z = (x, t, s) can be calculated 

G(z,z')=j: Mz) t x{Z ' ) • (38) 

In momentum space G is found to consist of two parts, G = G c hi ra i + G mass i ve . For G c hi ra i one 
finds 

<W«, cc / ^J^l e ^-^-^ . (39 ) 

This is the 2-dimensional massless propagator of a chiral fermion. The prefactor which is left 
out in (|39|) gives the s-dependence of the propagator. It reproduces the result we found earlier 
that it is bound to the wall with an exponential fall off when going away. A similar form for 



the propagator exhibiting its chiral structure was obtained in [^, 15 1. In addition to the chiral 
propagator also the one for the heavy modes can be computed 

G m assive(Z} % ) / TZ ' ^ ' A* 1; 2, 3 (40) 



with 



M=| mi2[S) I (41) 

1 -m(s n 1 



2 This is even more so for the following reason. In ]4^] it was pointed out that on a domain wall one should 
find the covariant anomaly. This is indicated by the strength of divergence of the Chern-Simons current. On the 
other hand, the anomaly of a 2-dimensional chiral theory is the consistent one as required by the Wess-Zumino 



consistency condition 50 . In 51 it was shown that by adding an extra term the consistent anomaly can be 
made covariant. It may be expected fl9|] that such an extra piece can be provided in the 2-dimensional system 
from the effective Chern-Simons action in 3-dimensions. 
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where 

k + ki 



2 2 m(s)m(s') + ik2(m(s') — m(s)) — uiq . (42) 
k 2 + Tn 

This describes a massive fermion. By analytic continuation to imaginary time the propagators 



m 12 



(|39|) and ( |40l) can be written in euclidean space where they will be denoted by S E and can 
such be used for the computation of the effective action in background gauge fields A p which 
in 1-loop is 

S eff [A] = l -J d 3 zd 3 z'A,(z)V^A u (z') (43) 

where 

= TrfS E (z, z')YS E (z } z') . (44) 



The expression in (|4"3"D can be calculated in the mo — » oo limit taking the gauge fields to be 
slowly varying. One finds the effective action to consist of two parts 

Seff = + Sffp 1 (45) 

with a Chern-Simons piece 

S?f f = J d 3 zsign(s)A fl d p A u (46) 

with fj,, v, p — 1, 2, 3 and a chiral piece 



S eff al = -^Ta / d 3 zd 3 z'A a (z)A b (z' f V + J V [1 - f(y, s - s')} 2 sech 2 (m s)sech 2 (m s') 

(47) 

where y = (x — x', t — t'), y a * = ie^yt is the dual of y a , the 2-dimensional index a = x,t and 

e — mo ||«|| 

f(y, s — s') — — — — [II^H cosh(m (s — s')) + (s — s') sinh(m (s — s'))] . (48) 

The Chern-Simons part S^ff is the known result . We see, however, that there appears 
an additional piece SfjJ al which has not been computed before. As we will see now this piece 
is the extra part as advocated in [49, 51 1 to find the covariant anomaly on the domain wall. 
We compute the currents by a gauge variation 

J » - If <«> 

J, = -^ign{s)e, pu d p A» + ±5, a e ab A b 5(s) . (50) 



From the current in (|50D we can see that we have obtained the Chern-Simons current as the 
first term. Its divergence comes from the different sign of the mass on the two sides of the wall. 
There is an additional piece in the current, residing exactly on the wall. It is this piece that 
gives an extra contribution to the divergence equations, making the anomaly on the wall the 
covariant one. 
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(it) We now proceed, discussing the second criticism of the Callan-Harvey computation, 
the regularization dependence of the Chern-Simons coefficient. For this purpose we describe a 
computation of the Chern-Simons current imposing a lattice regularization. We will work on 
an infinite lattice and keep the odd dimensions d = 3, although the results hold for arbitrary 
odd dimensions. We are interested in the low energy coefficient c of the Chern-Simons action 
Tqs which is obtained when the heavy fermions are integrated out. This leads to an effective 
action S eff = cT C s, 

Tcs = W / d^xA^dyAp . (51) 

The coefficient c is dimensionless and the Chern-Simons operator will therefore not decouple 
for large fermion masses. A discussion that heavy fermion masses may not really decouple from 



the low energy physics is given in [53|. The Chern-Simons coefficient c can be computed from 



the low energy portion of the vacuum polarization graph in fig. 4. One obtains 
i d f d^p 

C= 3! W &T J (2^)3 Tr [ S, ( p ) A A*(P'^~9)S , (p-?)A / ,(p-g,p)]| 9=0 . (52) 



Here S(p) is the free fermion lattice propagator which will be specified later and A M is the photon 
vertex. The integration in ( |52|) is understood to be taken over the 3-dimensional Brillouin zone. 
A first observation is that due to gauge invariance the photon vertex may be replaced in favour 
of the fermion propagator via the Ward identity 

A,(p,p) = -t^-S~ 1 (p). (53) 

OPv 

Upon differentiation with respect to d/dq p , the coefficient can be written as 



I {[StoW 1 ®] [S(p)duS- 1 ( P )} [S(p)d p S- 1 (p)}} . (54) 



The free lattice propagators S contains the lattice momenta sin(p). Therefore the lattice integral 
([54] ) appears to be quite horrible to compute. By exploiting its topological properties its 
calculation will become tractable, however. The topological significance of the above integral 
can be seen by noting that S 1-1 may be generically written as 

S~ x = a(p) + ib(p)a 

= N(p) [cos(|0(p)|) + i§asm(\9(p)\) 

= N(p)V(p) (55) 

where 



N(p) = V a 2 + b(p)b(p) , 0(p) = 6arctan(|6|/a) . (56) 

Quantities with an arrow denote a 3-dimensional vector and with a hat the corresponding unit 
vector. In this notation V(p) is seen to be a 2 x 2 unitary matrix. The integral ( |54"|) does 
not depend on N(p) provided that S" 1 ^) does not vanish. Thus S and S^ 1 may be replaced 
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S(p) 




S(p + q) 



Figure 4: The vacuum polarization graph. The low energy portion of this graph leads to the 
Chern-Simons coefficient to be calculated for the anomaly. 

everywhere by V and V\ The matrix V describes a mapping from the torus T 3 to the sphere 
S 3 . The integral ([54]) is then nothing else but the winding number of this map. Consequently 
it can only take integer values up to some normalization constant. The 3-dimensional example 
we have worked out here can be extended to arbitrary dimensions. 

To be specific we now take the usual fermion propagator for Wilson fermions for a constant 
mass m, 

3 3 

S^i'P) = ia v sin ^M + m + r H (! ~ cos Pm) • ( 57 ) 

fJ.= l fX=l 

The winding number can only change where S" 1 = 0. This happens only for momenta 
at the corners of the Brillouin zone and when the ratio m/r = 0,2,4,6 as can be seen from 
(|8|). For the ratio m/r — ► ±oo, V(p) — > ±1 and the integral vanishes. Therefore the Chern- 
Simons coefficient will be zero for m/r < and m/r > 6. Otherwise the integral (pl[) is 
piecewise constant and gives only contributions for momenta in an infinitesimal region around 
the Brillouin corners at the given values of m/r. To evaluate the integral we then only have 
to pick up these contributions and sum over them. Since the change of the winding number 
happens at the corners of the Brillouin zone we may expand the sin and cos function which 
leads us to evaluate 

dc . , . s u d f d 3 v m — 2rk 

dm 

„• 3 / Q \ 

(58) 

Our final task is then a trivial integration to get the coefficient as 




3 




3 (59) 
k J \m-2rk\ v ; 
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We want to remark at this point that the above derivation can be extended for arbitrary 
dimensions. The integral (|5^ ) describes then a map from the torus T d onto the sphere S d 
and the homotopy classes of these mappings are identified by integers. Therefore the whole 
calculation of the Chern-Simons coefficient follows very closely the above discussion. 

We now want to apply the result found above for the domain wall model we are interested 
in. We plot in fig.5 the integer part of the integral (|54| ) as a function of m/r (dotted line). 
The dashed line corresponds to the chiral zeromode spectrum as discussed in section (2.1). A 
value of +1 means one chiral fermion with positive chirality, —2 means two chiral fermions with 
negative chirality etc. Clearly the Chern-Simons coefficient follows exactly the behaviour of the 
change of the chiral zeromode spectrum with m/r. 

One peculiar feature of the Chern-Simons coefficient is that it is zero for negative m/r. 
Since in the domain wall model r is positive and m has different signs on the two sides of 
the wall, this means that the Chern-Simons current flows only on one side of the wall. This 
is to be contrasted with the continuum analysis which reveals that the current flows with 
equal strength but opposite signs on the two sides of the wall. Furthermore, if we take for 
example m — r — 1 the value of the Chern-Simons coefficient is c = —i/2n for m > which is 
exactly twice the continuum value. Thus we find that the strength of the Chern-Simons current 



depends on the regularization used as already emphasized in |47| . Of course, the divergence of 
the Chern-Simons current across the wall comes out the same in both calculations giving the 
correct strength of the anomaly. It is quite remarkable that the appearance of the anomaly 
on the lattice holds also for values of the domain wall mass at the order of the lattice cut-off. 
The fact that the Chern-Simons current flows only on one side of the wall motivates even more 
Shamir's approach to use free boundary conditions with the signs of the Wilson coupling and 
the mass chosen such that we will have current flow. 

The result for the 3-dimensional system discussed above can be generalized to arbitrary 
dimensions. In d — 2n + 1 dimensions the number of chiral zeromodes bound to the domain 
wall is for 2k < |^| < 2k + 2, < k < d - 1 



d- 1 

k 



(60) 



with the chirality of the modes to be (— l) fc sign(m). The corresponding Chern-Simons current 
is 



J% s (lattice) 
J^ 5 (continuum) 

with the continuum Chern-Simons current in d = 2n + 1 dimensions given as 




(61) 



i(-l) n 



^ ^4 7r y« n !| m | e ^ Ql -- a2 ™+ 1 ^ aia2 '''^ 1a2 ™ Q2 ™+ 1 • 
The above scenario can be implemented on a finite lattice of size L 2 L S . The euclidean lattice 
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Figure 5: The integer part of the Chern-Simons coefficient (dotted line) as a function of 
m/r. The dashed line indicates the chiral zeromode spectrum. +1 means 1 chiral zeromode 
with positive chirality, —2 is two chiral zeromodes with negative chirality. For m/r < and 
m/r > 2d the spectrum contains no chiral zeromode and the Chern-Simons current vanishes, 
accordingly. 



action for free domain wall fermions is given by 



Z 

J- r 



Z,fJ, 



(63) 



Z,fJ, 



with z = (x, t, s) and /i = 1, 3. We can couple external abelian gauge fields to the fermions in 
close analogy to the Hamiltonian ([H]) by making the finite lattice differences gauge covariant, 
see eqs. ( p9|j30D . Note that at this point the gauge fields are purely external. How to write 
down a dynamical gauge field action will be discussed in the next section. The lattice current 
is obtained by a gauge variation of the action and reads [TJ 



Jz 



+ 



r p 



We will choose a t-dependent external gauge field 



z,fi=2 



exp{— iq 



— £ cos(— t-1 

Z7T L 



} 



(64) 



(65) 



and make E 1 in order to stay in the low energy regime and below the critical momentum 
k c . The U's in the other directions have been set to one. 

As we are considering free fermions in an external gauge field background, the matrix 
elements of the current can be computed from the inverse fermion matrix using standard nu- 
merical techniques like Conjugate Gradient for the matrix inversion. Note, that in this way 
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the computation of the current is -up to rounding errors- exact. In particular, no simulation 
is involved. 

Since the width of the wavefunction in the extra dimension is finite, see fig.l, we will sum 
the current over the range in s corresponding to the support of the wavefunction. Then we 
evaluate the divergence 

< diji >= diji(t,x,s) ,i = (t,x) . (66) 

According to the discussion for the Hamiltonian formalism we expect the anomaly equation 
to be satisfied 

< 9di >= ±%-E eff {t) (67) 
where the effective electric field E e ff for small E is given by 

E eff = ^^E sm(t-l) (68) 

L 

and the sign is determined by the chirality of the mode. Choosing first the Wilson parameter 
r = and a charge of q = +1, it is found that the divergence of the current ( |6~6|) vanishes. This 
is perfectly consistent with the fact that for r = the doubler modes are still in the spectrum, 
cancelling the anomaly. 

Turning the Wilson parameter on should change the picture. We expect the doublers to 
become decoupled and the anomaly equation to be satisfied. The divergences of the currents 
computed for charges of q = 3, 4, 5 and chirality +1, +1, —1 are shown in fig. 6. Each of them 
follow the anomaly equation individually up to a few %. Taking the sum of them the total diver- 
gence vanishes which corresponds to the anomaly cancellation as predicted by the Pythagorean 
relation 3 2 + 4 2 — 5 2 = 0. In particular, the numerical results confirm that the Chern-Simons 
current flows only on side of the wall fl54 |. Thus we see both, the individual anomaly for a 



fermion with a given charge and the cancellation of the anomalies if the fermions are in an 
anomaly free representation. 



5 Coupling to gauge fields 

Despite all the positive results we found in the previous sections the crucial question remains, 
whether we can keep the chiral zeromodes decoupled from each other when the gauge fields 
are made dynamical. The danger is that the gauge fields might induce an interaction between 
the two domain walls. We will discuss two proposals of coupling gauge fields to domain wall 
fermions. In a first attempt, the gauge field is 5-dimensional. The gauge couplings are chosen 
differently for the 4-dimensional fields and in the fifth direction. The hope in this approach is 
that by tuning the 4-dimensional gauge coupling to zero and making the coupling in the fifth 
direction strong at the same time, physics is confined to 4 dimensions. One then might perform 
the usual continuum limit in 4 dimensions while keeping the chiral structure of the theory. 
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Figure 6: The divergence of the (l+l)-dimensional gauge current in an external gauge field. 
The system consists of fermions with charge q = 3, 4, 5 and chirality +, +, — , respectively. Each 
of the currents obey the lattice anomaly equation individually. The sum of the currents vanish 
due to the anomaly cancellation. 

We will see, however, that this path does not lead to the desired result, namely a 4- 
dimensional chiral gauge theory. The failure of this approach suggests that one has to keep the 
gauge fields strictly 4-dimensional from the very beginning. In order not to couple both walls, 
the gauge fields have to be switched off on one of the domain walls such that only one wall is 
gauged. We will thus have a situation with identical 4-dimensional gauge fields on a number of 
s-slices around a wall and U — 1 in the complementary region around the anti-wall. Clearly, 
at the boundary where both regions meet we loose gauge invariance. This can be repaired by 
introducing a scalar (or Stuckelberg) field at this boundary. However, as we will see, there exist 
light mirror fermions living on the boundary. They can interact with the chiral zeromode on 
the gauged domain wall thus rendering the theory vectorlike again. 

5.1 Dimensional reduction 

In this first way to couple the gauge fields the gauge interaction is split into a purely 4- 
dimensional part and a piece that contains the gauge fields in the fifth direction. Both parts 
are equipped with different gauge couplings. For the possible phase structure and the question 
of whether one may obtain a chiral gauge theory we will consider a 5-dimensional system. The 
reason for making a departure from our 3-dimensional setup will become clear later. For the 
discussion it is sufficient in a first step to consider the pure gauge theory alone. We will choose 
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U(l)-gauge fields. The action becomes 

s(u) = E ^E MK^+^U^) + & E ^(c/^kv^TK") ] • (69) 

Here /?4 = l/#f is the 4-dimensional and fa = l/gf the 5- dimensional inverse gauge coupling 
squared. We have denoted with U the gauge fields living in 4 dimensions and with V the ones 
in the extra dimension. Thinking of s more as a flavour space than an extra dimension, we have 
denoted the s-dependence with a superscript on the fields. Note that the role of the l^-fields 
look very much as that of scalar fields coupling to 4-dimensional gauge fields in some peculiar 
flavour space. We are interested in the properties of this 5-dimensional gauge theory in the 
limit of large f3± corresponding to the continuum limit in 4 dimensions. 

Letting the gauge fields be U X;/X = e l9x ^ and V = e %e ' x ^ the path integral is 

z = J ve x ^ve' X41 e- s ^ e '^ ] . (70) 

We consider this model in the meanfield approximation [55, |56| which is obtained by inserting 

1 = J T>v XjfM [ Va exp J ^ [a XyjJi (v XjfM - 9 X ^)\ \ (71) 

in the path integral. In ([71]) the mean field v is a complex variable living on the links of the 
lattice and a an auxiliary field. There is an analogous expression for 9' . Now the integration 
over the original link variables 9 and 9' decouple and the path integral becomes 



f Vv f Va exp J ]T ( ft ]T Reu^ + ft ]T Rev x ^ 5 j + W(a) + W(a' 



) — v a — v'a' 

(72) 
where 

W(a)= I** d9 x ^ a6 ^ (73) 
Jo 

and v x/M/ denotes the product of the v's around an elementary plaquette in 4 dimensions and 
v x ,fi5 the corresponding expression in the fifth direction. A saddle point of the action in ([72]) is 
obtained by solving 

SS _ §S_ l 

dv " ' dv' " Cja\ 

dW = y . dW = v i y'^> 

da ' da' 

Searching for translational invariant solutions by setting v and v' to constants, the explicit 



equations at tree level are [56 



a = 2(d - 2)f3 4 v 3 + 2(3 5 v' 2 v + 2(3 4 v 

a = ft (2(d- l)v 2 v' + 2v) (75) 



and 

V = I (a) ' " - Jo(a') 



v = m (^) 
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Figure 7: The phase diagram as obtained from a numerical simulation of the 5-dimensional 
U(l) gauge theory for an 8 4 lattice. Only the portion for large (3^ is relevant for the continuum 
limit. 

where 

I k ( a ) = - r d6cos k (6)e acosW . (77) 

7T JO 

In principle there appear three different phases see f56[ and fig. 7. However, we are only 
interested in the large f}± limit and will not discuss the phase with v = and v' — 0. For large 
/?4 one finds a layered phase with v > 0, v ' = and separated from it by a phase transition at (3% 
a symmetry broken phase with v > and v' > 0. The layered phase is characterized by the fact 
that a charged particle will only move along the layer and can not hop between layers. Therefore 
if one would live inside a layer one would basically only experience 4-dimensional physics. In 



56| it was argued that such a phase only appears when the lower dimension is greater than 
two. Thus an investigation of a 3-dimensional system would simply miss the layered phase and 
might lead to wrong conclusions. In fig. 7 we show the phase diagram of the 5-dimensional pure 
U(l) gauge theory. It is obtained by means of numerical simulations on an 8 4 lattice. As is 
predicted by the meanfield computation, there are indeed two phases at large f}±. One of them 
has v > ans v' > and corresponds to the symmetry broken phase. For /3 5 < /3| there appears 
the layered phase which is the promising phase to obtain 4-dimensional physics. Taking the 
1-loop corrections to the mean field equations into account, the phase diagram remains stable. 
The fermions are included via the action 

X,S fj,= l 
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+ \ Et£ [V? WC 1 + 75) + V x * s 5 s ^ & , s (l - 75) - 25s, s] r x (78) 

where we have separated the 4- dimensional part from the part in the fifth direction and set 
r = 1. Again we have written the s- label as a "flavour" label and ^5 denotes a unit vector in the 
5-direction. In EJ| it was shown that the fermions lead only to a slight modification of the mean 



field equations and that the phase structure depicted in fig. 7 remains stable when fermions are 
taken into account. The interesting question, of course, concerns the dynamics of the fermions 



in the layered phase. One might get a hint by an inspection of the action itself EH]. In the 



layered phase v' = corresponding to set V in the action ( ff8| ) to zero. But then the action 



looks like a normal 4-dimensional Wilson action in each layer. Indeed, in [24| the fluctuations 
to the tree level saddle point solution were computed to get the fermion propagator. As a result 
one finds for every s-slice the 4-dimensional Wilson propagator. This means that the theory 
-though truly 4-dimensional- is vectorlike in each s-layer. Therefore we have to throw away 
the region /3 4 3> 1 and (3$ < (3$ as a possible corner of the phase diagram for a construction of 
a chiral gauge theory from domain wall fermions. 



In |27fl it could be proven by the hopping parameter expansion that for /J5CI the fermion 
propagator is parity invariant, with the parity transformation defined such that for Xj — > —Xi 
with x being the 4-dimensional coordinate 

V> 70^ , ^ -> tPlo ■ (79) 

Therefore there is a complete symmetry between left and right handed modes. Thus the theory 
is completely vectorlike in each layer strengthening the above conclusions. 

What remains is the possibility that (3$ > {3$ |2"5fl , i.e. v > 0, v' > 0. In this situation we 



have symmetry breaking. In each s slice there is a G — U(l) (g> U(l) symmetry such that the 
total symmetry group is G La when the extent of the system in the fifth direction is L s . For 
v' > this symmetry is broken to its diagonal subgroup G. Then we are left with L s — 1 massive 
gauge bosons with a mass mc oc v '. However, there will remain one massless gaugefield which 
does not depend on s. This gauge field couples equally to the modes at the domain and the 
anti-domain wall. Therefore, although we still would have the zeromodes on the domain walls, 
they can communicate via the massless gauge boson. We expect therefore the model again to 
be vectorlike. 



5.2 Waveguide model 

The previous section showed that starting with 5-dimensional gauge fields does -most probably- 
not lead to a chiral theory via dimensional reduction. Obviously the gauge fields have to be 
strictly 4-dimensional from the very beginning. This becomes a more natural point of view if 
one considers the extra dimension as a flavour space with a somewhat unusual flavour matrix. 
On the finite lattice with two domain walls not all s-slices (flavours) can be gauged. This would 
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immediately lead to the possibility that the zeromodes on both domain walls communicate, 
rendering the theory vectorlike. 

One might therefore try a scenario where only one of the domain walls is gauged and 
can interact with the gauge fields. This corresponds to taking a number of s-slices around a 
domain wall and put identical 4-dimensional gauge fields on each of these slices leaving the 
complementary s-slices ungauged. That the gauge fields have to be identical is again motivated 
from the flavour space picture. In this way both walls are completely shielded from each other. 
However, following the above prescription one notices that at the s-layer where a gauged s-slice 
meets an ungauged layer gauge invariance is broken. 



One possibility is that one does not worry about broken gauge invariance [p7[l . Indeed, 



several chiral fermion proposals exist that start with broken gauge invariance JT^. It is then 
hoped that gauge invariance is restored in the continuum limit. However, we want to follow 
a path where gauge invariance is kept. Its restoration is actually very easy. One just has to 
remember how the Yukawa-coupling in the Standard Model is made gauge invariant, namely 
by introducing scalar (or Stiickelberg) fields. From this point of view, the model with broken 
gauge invariance can be thought of as the gauge invariant model in the unitary gauge. Thus, 
both descriptions are actually completely equivalent M. The proposal to couple gauge fields in 
the domain wall model is 

• keep the gauge fields strictly 4-dimensional, 

• gauge only a number of s-slices around one domain wall, 

• introduce scalar fields at the boundary of the gauged region. 

The gauge fields are then confined to a region in s between the boundaries where the scalar 
fields live. Thinking of electromagnetism this resembles a situation where the electromagnetic 
waves are trapped in a waveguide which suggests the name "waveguide model" that we will use 
from now on. 

The introduction of the scalar fields gives rise to a coupling between the fermions and these 
fields. It is then natural to equip the resulting interaction with a Yukawa-coupling y. To get 
a feeling for the physics of the model proposed, one could first switch off the Yukawa coupling 
and set the gauge fields to one. In this situation the zeromode spectrum can be obtained again 
from the Hamilton formalism. The result for a 3-dimensional system is shown in fig. 8a. 

The situation resembles a combination of the domain wall model and the boundary fermion 
model. Setting y = cuts the waveguide region completely off and produces open boundary 
conditions. Due to our discussion in section (2.2) we expect to find chiral zeromodes where the 
mass term is negative. This is indeed what happens. We find the original chiral zeromode on 
the domain wall ipi and a chiral surface mode -denoted in the following as "mirror fermion" - 
with opposite chirality xr a t the boundary. A complementary picture is obtained in the s-region 
outside the waveguide. 
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Figure 8: The zeromode spectrum in the waveguide model. The boundary where the scalar 
fields live are indicated by the shaded stripes. The position of the domain walls are given by 
the vertical bar. Fig.8a corresponds to y = and fig. 8b to y > 0. 
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Staying in a meanfield setup by giving the scalar field a vacuum expectation value v, we 
show in fig.8b the spectrum for a non-zero value of y. Now the mirror mode can leak out of the 
waveguide region and can combine with the mirror mode from outside the waveguide region to 
form a Dirac particle. The mass of this particle is expected to follow the usual perturbative 
relation rrif = yv . 

Let us make clear what we would like to achieve with the waveguide model. Clearly one 
wants the chiral zeromode on the domain wall to survive. At the same time, the mirror modes 
at the waveguide boundary should be heavy and therefore decouple from the spectrum. This 
scenario looks hopeless if the mirror fermion mass follows the perturbative relation. Approach- 
ing the phase transition to the symmetric phase, v (in lattice units) approaches zero. Then 
also the mirror fermion becomes massless and can consequently not be decoupled. Moreover, 
in the symmetric phase the mirror fermions are expected to be massless since there the vacuum 
expectation value vanishes identically. We would therefore end up with a variant of Montvay's 
mirror fermion model []5S] . 

What gives hope for the waveguide model is its resemblance to Yukawa-models on the 
lattice. In fact, concentrating only on the waveguide boundary one might consider the model 
as being just of the Yukawa-type. Such models have been investigated intensively in the last 
years. The most surprising and unexpected outcome of these investigations was the discovery 
of a new phase at large, non-perturbative values of the Yukawa-coupling |32|, |33| . 

When the Yukawa-coupling is very large, the fermions can easily combine with the scalar 
fields and form massive bound states even in the symmetric phase where the vacuum expectation 
value is zero. The masses of these states can be at the order of the cutoff such that they decouple 
in the continuum limit. What made the Yukawa models not successful for a regularization of a 
chiral gauge theory was that the continuum limit in the strong coupling region did not resemble 
the Standard Model at all. The fermion spectrum in the continuum limit consisted of only a 
neutral Dirac fermion. The necessary charged fermion that would couple to the left handed 
gauge field is simply missing and does not appear in the spectrum. What is important is, that 
on the other hand it was easy to give the neutral fermion a large mass. 

For the domain wall model the situation is different in that we still would have the chiral 
zeromode on the domain wall. Here the strong coupling phase would just serve to make the 
modes at the waveguide boundary heavy. The real physics -and hopefully a chiral theory- 
would appear along the domain wall. The main question in the domain wall approach to chiral 
fermions on the lattice is then whether in the waveguide realization a strong coupling phase 
exists. 

This makes it necessary to explore the phase diagram of the waveguide model. In doing so 
we will impose some simplifications. First, we set the gauge fields U = 1. This is certainly 
justified as in the continuum limit the gauge coupling becomes small and the gauge fields might 
be treated perturbatively. Furthermore, the investigations of Yukawa models on the lattice 
revealed the strong coupling region leaving out the gauge fields, too. Therefore it appears 
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Figure 9: The phase diagram which would make the waveguide model successful. 

to be sufficient to study the system with fermions and scalar fields alone. Secondly, we will 
go back again to a 3-dimensional system for the numerical simulations that will be presented 
below. Experience from earlier work suggests that they resemble 4-dimensional systems quite 
closely p9] |. In particular the strong coupling behaviour could also be identified in this lower 



dimensional models. Thus the appearance or non-appearance of a strong coupling region in 
the 3-dimensional model would strongly point towards that such a region may or may not exist 
also in higher dimensional models. 

A sketch of the phase diagram which would make the waveguide model successful is shown 
in fig.9. There are two Yukawa-couplings. First there is the coupling we introduced at the 
waveguide boundary to couple the scalar fields directly to the fermions. Secondly, due to the 
exponential fall off, the domain wall fermion will have an -exponentially suppressed- overlap 
with the modes at the waveguide boundary and hence couples to the scalar field with a strength 
2/exp (— m L s /4:). (The distance of the waveguide boundary to the domain wall is L s /4, see 
fig. 8). The phase diagram in these two couplings may have the following phases. 

• FM In the ferromagnetic phase the symmetry is spontaneously broken and the scalar 
field assumes a vacuum expectation value v > 0. If the model would be successful, there 
should be two regions in the FM phase. One with weak coupling behaviour mp pa yv and 
a strong coupling region with heavy fermions. 

• PMW This is the weak paramagnetic phase. Here the symmetry is restored and the 
vacuum expectation value v — 0. In this phase the mirror fermions living at the waveguide 
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boundary are massless and do not decouple from the spectrum. Therefore it is not possible 
to construct a chiral gauge theory in this phase. 

• PMSi A strong paramagnetic phase. This is the phase we are looking for. Here due to the 
large Yukawa-coupling the binding of the fermions to the scalar fields at the waveguide 
boundary becomes so strong that they combine and form massive bound states. At the 
same time the chiral zeromode on the domain wall remains massless. Therefore the mirror 
fermions are expected to decouple and we are -hopefully- left with a chiral gauge theory 
on the domain wall. 

• PMS2 Again a strong paramagnetic phase. However, this time the extension of the system 
in the extra dimension is so short that the overlap of the chiral zeromode on the domain 
wall with the waveguide boundary becomes non-negligible. Therefore the domain wall 
modes can also combine with the scalar fields getting massive in this way. Clearly also in 
this phase no chiral fermions can be obtained. 

In the following, we will use PMW to indicate a weak coupling behaviour of the fermion mass, 
to/ oc v. PMS will indicate strong coupling behaviour with massive bound states. Note, 
that in this way PMS and PMW will not necessarily stand for weak and strong values of the 
Yukawa-coupling y. 



5.3 The action of the waveguide model 

After the qualitative discussion of the previous subsection, we want to become more concrete 
and will construct the lattice action suitable for eventual simulations. As we discussed above, 
we take the same gauge field on all s-slices inside the waveguide. The index for the extra 
dimension will be chosen as a superscript to show the resemblance to the flavour space picture. 
We define gauge transformations on the fermion field as 



* s x -> g x V s x , * x -> * x9 l seWG, 



WG = {s : s < s < s' } (81) 

with g x in a gauge group G. The detailed choice of the boundaries so and s' is not very impor- 
tant, provided they are sufficiently far from the domain wall that the zeromode is exponentially 
small at the waveguide boundary. For symmetry reasons we shall choose s = (L s + 2)/4 + 1 
and s' = (3L S + 2)/4, such that the right handed mode at s = L s /2 + 1 is located at the center 
of the waveguide, see fig. 8. With this choice we have to take L s — 2 a multiple of four. 

Having made this division into a waveguide and its exterior, we note that the model has a 
global G x G symmetry: 

* x -> 9 V X , X -> WJ, seWG, 
V x -> h*%, K -> Kh\ s^WG. 
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With our choice for the position of the waveguide boundary, there is a symmetry involving 
parity plus a reflection in the s-direction with respect to the plane s = so — \ — L s /A + 1, 

K - 7^pf +2 ~ S , (83) 

with Px = (—Xi, • • • , —Xd-i,Xd) the parity transform of x. 

The layers, where gauge invariance is broken are now sq — 1 to sq and from s' to s' + 1. 
As discussed in the previous section, this will be repaired by the introduction of scalar fields 
(or Stiickelberg field) V at the boundary of the waveguide. Alternatively the scalar field V 
might be interpreted as a component of the gauge field in the extra dimension. In this sense 
the waveguide model is an improved variation of the actions flB5j , |7j|) where the gauge fields have 
been made 5-dimensional. It corresponds to making the gauge coupling in the fifth dimension 
/?5 = oo except at two s-slices, namely at s = sq and at s = s' . We obtain the gauge invariant 
action (suppressing the even dimensional index x) 



seWG s(£WG 

- y(F°"VPi* S0 + FV^ijf 80 " 1 ) - y(W'°V^P L ^ s ' 0+1 + W' 0+1 VP R y s '°), 

where we have supplied the Yukawa term with a coupling constant y as anticipated. Note that 
we take the same scalar field at both waveguide boundaries. Since we have chosen r = 1 we 
have written projectors in the hopping terms in s, Pr(l) = 1(1 + ( — )7s)- $(U) = l^,d x {U) is 



the usual gauge covariant lattice Dirac operator , see ( p9| , pOD . The field V x G G is the scalar 



field, which can be thought of as a (radially frozen) Higgs field, and which transforms as 

V x - hV x gl (85) 

The transformation given in eq. ([83]) remains a symmetry if V transforms as 

V x - V£ x . (86) 

Having introduced the scalar field, it is very natural to add a kinetic and an interaction 
term for it. It is expected that through renormalization these terms would be generated auto- 
matically. Since we have the scalar field to be radially frozen, we choose the standard lattice 
form of the non-linear a model 

Sv = - K Y Tr(V x U XtlM V% +lt + h.c). (87) 

X,fl 

The scalar hopping parameter k is proportional to the bare mass squared. For the pure scalar 
theory we will have a broken phase for k > k c as well as a symmetric phase for k < k c . 
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5.4 Mirror fermion representation of the model 

To get more insight into the mode spectrum of the model and to make the interpretation of the 
surface modes xl,r m fig-8 as mirror fermions more plausible, we rewrite the action as follows. 
Relabel the right and left handed fermion fields, ^f s RL = Pr : l^ s as 

Xl = xk = *S"*, (88) 

and the same for $?R t L = ^Pl,r (note the reversal of L and R). The new label t runs from 1 to 
L t = L s /2. In fig. 8 we have indicated this new labeling for the zeromode wave functions shown 
there. With our choice for sq, s' and L s we can define a domain wall mass for both fields ip 
and Xi which is a step function in t satisfying, 

jjt = m *o-i+t = m s -t (gg) 

With this relabeling the two domain wall zeromodes will reside in the Dirac fermion field if), 
whereas the waveguide boundary zeromodes will reside in x- After substituting eq. (|88|) into 



eq. fl8"4|) with U = 1, the action turns into 



Lt 



t=i 

L t -1 



t=l 



- ytiVPL + V^P^x 1 -y^ L \V^P L + VP R )^ L \ (90) 

In this form, the action resembles that of an L r flavour mirror fermion model in the fashion 
of ref. p8"| , with ip the fermion and x the mirror fermion field. In fact, for L s = 2 the hopping 
terms in t are absent, /Z* = and the model reduces to the mirror fermion model of ref. |58| with 
equal Yukawa couplings for the fermion and the mirror fermion, and a vanishing single-site mass 
term. For L s > 2 the model has a more complicated mass matrix (i.e. non-diagonal couplings 
among the flavours s or t) and if the model is going to be more successful in decoupling the 
mirror fermion than the traditional mirror fermion approach, it must come from this mass term. 

The mass matrix for the L t flavours in the model is not diagonal but this can be remedied 
by more rewriting. First we expand the fermion fields in a plane wave basis, which diagonalizes 
the Dirac operator and Wilson term, ip s x = J2 P e lxp ipp, tp S x = J2 P e~ lxp tjj s p . Here J2 P 1S a normalized 
sum over the momenta on the d dimensional lattice, Y^ P 1 = 1. Then we can write, 

Lt _ _ 

s i>x = EEO^P^ + ^pXp+^K+^^ + ^K+^xJ 

Lt-l 



E + 4 + VJ) + E &Ui + m) 

t=l t 

yEfe^^ + ^^^xS+^^Pi + v^vf), (9i) 



pq 
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with ftp = YJ 7^ sin(p M ) , w p the diagonal form of the Wilson term, w p = — cosQo^)) and 

V p the Fourier transform of V x . For y = the action has the schematic form 

*-<«U :/)(:)• 

with M a (p dependent) matrix in flavour space, which can be read off from eq. fl9~l~l). This 
action can be diagonalized by making unitary transformations on ip and \i 



t f = G) tX \ I s = tG t f, 



(93) 



such that G^ s M st F tg = fifSf g . The matrices F and G are eigenf unctions of M^M and MM^ 
respectively, labeled by the index /: 

(M^M) st F tf = \fi f \ 2 F sf , (MM^) st G tf = \fi f \ 2 G sf . (94) 

For suitable choices of the phases of the eigenfunctions, the fi/s are real. Substituting the 



mode expansion (|93T) into the action ( |DTD with the momentum label restored, we arrive at 



f=i p 

-V E %G%{V p _ q P L + V^ p P R )Gl g e g +uJ f p Ff Lt (vl P P L + V p . q P R ) F^ . 
fg,pg 

In this representation of the model, it is seen that all fermion modes uo? and interact with 
the scalar field, but that their effective Yukawa- coupling is determined by the magnitude of 
their wave function at the waveguide boundaries t — 1 and t = L t . This leads then naturally 
to the introduction of the two Yukawa-couplings we used in fig. 9 to parametrize the phase 
diagram. There we use the value of the wavefunction of the mirrormode at the waveguide 
boundary which gives just y and the value of the domain wall fermion at the boundary which is 
y e -m-oW4 p or y — q tn e model is seen to describe free, degenerate fermions and mirror fermions 
with momentum dependent mass \J p (for \rf ^ 0, the eigenstates are + £^ and ujj — £f). 
Exactly one flavour, which we denote with / = 0, has //° = (up to terms exponentially 
suppressed in L s ) for \p\ < p c , where p c is the critical momentum defined in sect. 2.1. For r = 1 
and m close to 1, the critical momentum is p c ~ y/2. All other p,f and also p, p for p outside 
the critical momentum region, are 0(1) in lattice units. 

This shows that for y = and momenta \p\ ~ p c , the model contains a massless fermion, 
u°, as well as a massless mirror fermion, £°. All other modes (/ 7^ 0) as well as the species 
doublers have a mass of the order of the cutoff. The species doublers of the zeromode / = 
are massive because n° p is 0(1) for momenta with p^ = ±7r. 

As was discussed already in section 5.2, fig. 8 shows the t-dependence of the zeromodes F t0 
and G t0 of the fermion (indicated by ip in the figure) and mirror fermion (indicated by x) f° r 
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the smallest momenta \p\ = it / L <C p c . It shows that the zeromode for the fermion is sharply 
peaked at t = (L t + l)/2, i.e. at the domain wall and the zeromode for the mirror fermion is 
localized at the boundary, at t = L t . The non-zeromodes, which are not shown in this figure, 
are not localized. 



5.5 Large y expansion 

The previous sections introduced the action of the waveguide model and through a mode 
expansion exhibited its mode spectrum. In addition to the chiral zeromode on the domain wall, 
new mirror fermion modes appeared at the wave guide boundary. It can be expected that the 
mirror zero modes found for y — 0, which are depicted in fig. 8a, remain light and proportional to 
the vacuum expectation value v for small but non-zero y. The interesting question is whether 
for large values of the Yukawa-coupling the zeromodes at the waveguide boundary become 
massive and eventually decouple. 

As fig. 8 and the discussion of the previous section suggests, things may be simplified by 
regarding a system with one waveguide only and to consider the effective scalar-fermion model 
at the waveguide boundary alone. To get some insight into the spectrum at large y, one can 
try to perform a 1/y expansion of the so truncated model [B0|. Again the gauge fields will be 
set to U = 1. The 1/y expansion is obtained by first relabeling the fermion fields in s, while 
suppressing the even dimensional index 



(96) 



By rescaling the field ip 

S F = J2 (Xa@Xs + Xs a ( S )(~ A + m ~ !)Xs) 
s 

+ Yl (xlsXRs+i + Xr s +iXls) + ^ L y^R + ^PrVHl 

s 

+ y/a A +m- l)x s ' o + + m - + mp^ip, (97) 

where ip R>L = P RiL if> , xr,l = Pr,lX, 

a = — and a(s) — 1 — S a j . (98) 

y s '*0 

The sums over s are to be taken in a region around s' neglecting all other defects in the model 
which are thought of to be far away. In particular the domain wall mass is constant. 

In this form it is seen that for y — ► oo, the ip and the \ fields decouple. For the x-fields at 
the waveguide boundary one may derive the Dirac equation and look for solutions that simul- 
taneously satisfy the massless Dirac equation of the lower dimensional system. Imposing plane 
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waves in the even directions the analysis proceeds very similar to the discussion in section 2.1. 
The Dirac equation reads in momentum space 



#(p)x. + a(s)(m - 1 - F(p)) X s + PrXs+i + PlXs~i = 0. (99) 

Solving eq . (|99"|) for the x-fields, one arrives at normalizability conditions that dictate the chiral 
zeromode spectrum on the waveguide boundary |BD| . It turns out that in 4 + 1 dimensions one 



finds 7 righthanded and 8 lefthanded zeromodes |Q inside the waveguide, which are bound to 
the waveguide surface. If the gauge fields are switched on, this leads to a vectorlike spectrum 
if one takes the zeromode at the domain wall into account. 

An important question for a further exploration of the phase diagram is whether the above 
result is stable against perturbations in a. To answer this question we will now compute the 
X-field propagator at the waveguide boundary to order a. One may start by integrating out 
the -^-fields in ( p7|) which leads to an effective action 



S x = (Xs$Xs + Xs a ( s ) (- A + m- l)xs) + E {X s PrXs+i + X s +i p lXs 

s s 

oo 

- «E X s ;(-A + m - 1) J2(-<[( V ^ p R + VP L W} n (V^P R + VP L )(A + m- l) Xa > 

n=0 

= Y.X s SfyXs> - «EtX/ (- A + m - WMPr + V X P L )[(A +m- l) Xs ' Q l 

s,s' £ 

+ 0{a 2 ). (100) 

Here is the free x-held propagator which can be computed exactly, see e.g. section 8. To 
order a, replacing V by its expectation value vl, the Dirac equation reads now 

iKP)Xs + a{s)(m - 1 - F(p)) Xs + PrXs+i + PlXs-i - av(m - 1 - F(p)) 2 8 s ^ Xs > o = 0. (101) 

For the symmetric phase with v = this is exactly the same equation as we found for y = oo. 
Therefore the result from the large-y expansion is, that in the symmetric phase we will find the 
same rich zeromode spectrum with equal numbers of left and right handed fermions as we had 
at y = oo. This is quite in contrast to the anticipated phase diagram, fig. 9, where we expected 
actually a symmetric phase with strong coupling behaviour. 

For v > which corresponds to the broken phase we have an extra piece in the Dirac 
equation at s' , av(l — m + F(p)) 2 . This gives a mass to the mirror modes of size v/y. We see 
that the weak coupling relation rrif oc yv is replaced by rrif oc v/y, i.e. y — > 1/y. In particular, 
the fermion masses are still proportional to v and remain light when approaching the phase 
transition, v — > 0. Passing the phase transition we reach a symmetric phase with a vectorlike 
zeromode spectrum. In the discussion of the large-?/ expansions we experienced a first surprise. 
Since in the large-?/ expansion all defects are thought to be very far apart from each other, we 
effectively work in the limit L s 1. This means that we are on the jy e - m oW 4 ~ Q-axis in the 
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phase diagram of fig. 9. There we hoped to have the PMSi phase at y 3> 1. Instead we find 
a second symmetric phase with weak coupling behaviour and an even richer -but vectorlike- 
zeromode spectrum than the one at weak Yukawa- couplings. 



5.6 Reduced model 

The mode expansion discussed in section 5.4 revealed more clearly than the original action ([84]) 
what the important modes of the waveguide model are. There is, of course, the domain wall 
fermion living on the gauged wall. Furthermore, there are mirror fermion like states that live on 
the waveguide boundary. These states can easily get light and remain in the spectrum. We saw 
that this can happen for weak coupling as well for very large coupling, which led to a revision 
of the phase diagram in fig. 9. However, at intermediate values of y, there might still be still our 
candidate phase, PMSi where we can hope to construct a chiral gauge theory on the domain 
wall at the end. It became clear that it is very natural to think of the model as a mirror fermion 
model along the lines of [|58| and that it would be the peculiar mass matrix which would be 
responsible for decoupling the mirror fermions in contrast to the normal approach |58] . 



As we will see below, numerical simulations in the waveguide model even in the quenched 
approximation become very demanding. On the other hand we could identify the important 
building blocks of the domain wall model, the behaviour of which make the model succeed or 
let it fail, uniquely. These building blocks are the domain wall mode and the waveguide mode 
inside the waveguide. Furthermore, these modes appear only in the low energy regime below 
the critical momentum p c . Above p c they become heavy and delocalize. 

One may therefore attempt to construct a reduced model that shows the characteristic prop- 
erties of the waveguide model with these ingredients. As we saw in section 5.5, the anticipated 
fermion spectrum for the reduced model is not realized at very strong values of y. Therefore 
the reduced model only makes sense under the assumption that at intermediate ^/-values the 
PMSi phase really exist. We will use for this "reduced" model the domain wall zeromode uj°, 
the mirror mode £° at the waveguide boundary, the scalar fields^ and the critical momentum 
p c . This leads to a model 

sr* = Y,H^° P + ^pQ + y E f P {vl P p R +v^ q p L )e q . (102) 

\p\<Pc \p\,\q\<Pc 

In this reduced model the domain wall fermion u° is completely decoupled from the mirror 
fermion £° and the scalar field V. The mirror mode is then coupled to the scalar field with 
a strength of the Yukawa-coupling y. We now have exactly the situation which might lead to 
a successful description of a chiral theory if the waveguide model would work: An unaffected 
chiral zeromode on the domain wall and a Yukawa-like theory on the waveguide boundary. The 
only question is, whether this Yukawa-like theory possesses the strong symmetric phase or not. 
One may argue that the above model is an oversimplification of the full theory. But one can 



give plausible arguments that this kind of reduction makes sense |25[ . We will see further from 
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the numerical simulations that the reduced model behaves almost the same as the full model. 
This is so at least at the values of y where the simulations discussed in the next section have 
been performed. We know, of course, from the large-?/ expansion, that the reduced model does 
not represent the full model at y 3> 1. Third, if in the reduced model no PMSi phase exists, it 
will also not exist in the full model. The only addition from the full model would be that the 
domain wall zeromode and the mirror mode start to interact. But then they either combine 
to form a massive Dirac particle or the effective interaction is so weak that we are back in the 
situation of the reduced model. The additional heavy modes in the model, in particular the 
modes above the critical momentum p c , do not play an important role as their coupling can be 
expected to be inversely proportional to their mass and hence is weak. 



6 Numerical simulations in the waveguide model 

Up to now we discussed the waveguide model at the two boundaries of the phase diagram for 
small and large Yukawa-couplings. In both these limits it seems not to be possible to construct 
a chiral gauge theory. This may, however, still be possible in the middle of the phase diagram 
at intermediate values of the Yukawa-coupling. In this region both, the weak coupling as well as 
the strong coupling expansion do not work. Therefore, in order to answer the question whether 
the PMSi phase exists one has to rely on numerical simulations. 

As already mentioned earlier, Yukawa models in 2 and 4 dimensions have been subject to 
extensive investigations in the last few years. In these studies it became clear that they could 
either show a PMSi phase [j5j| or not H6l| , depending on the choice of the Yukawa-coupling. 



Two examples of these kind of Yukawa models are the local coupling 

Sic = y y E^x(v x PR + v x *p L )ij x , (103) 

X 

and the hypercubic coupling 

She = y £ ~ Y,^M-i>Pr + V X %P L )^ (104) 

x 4 b 

where we only show the Yukawa-part of the action. In the case of the hypercubical coupling 
the sum over b means a sum over all points of a lattice hypercube surrounding the point x. The 
model with the local coupling is known to have a PMS phase whereas the one with hypercubical 
coupling does not. 

A form of the Yukawa-coupling with a momentum cut-off as in the waveguide model has so 
far not been investigated. Therefore the question whether a PMSi phase exists or not remains 
undecided. Unfortunately the special form of the Yukawa-coupling as given in the wave guide 
model and in its reduced form prohibit an analytical treatment for intermediate values of y. 
Obviously a numerical study of the 5-dimensional system would be very demanding. It turns 
out, however, that the strong coupling behaviour we are after already appears in 2-dimensional 
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Yukawa models. Even more, the phase diagrams of 2 and 4-dimensional models are very similar 
so that it is a reasonable strategy to search for the PMSi phase in the waveguide model in 2 
dimensions. 

This path has been followed in fl2"5| . Most results are obtained in the quenched approxi- 



mation, neglecting all fermion loops. This seems to be sufficient for the question whether a 
PMS phase exists. Although in the quenched approximation the phase transition lines shown 
in fig. 9 do not really exist, the model should show typical weak or strong coupling behaviour. 
The experience from earlier studies of Yukawa models suggests two approaches to identify or 
exclude a strong coupling symmetric phase. The first is the fermion mass spectrum and the 
second the eigenvalue spectrum. Both show characteristic behaviour if a strong coupling phase 
exists. The results of this numerical investigation will be discussed in the next two subsections. 

The (2 + 1) -dimensional model is simulated using U(l) scalar fields. The pure bosonic 
sector of the model consists therefore of a YY-model on the boundary of the waveguide region. 
One might worry that the XY-model we are going to look at is too different from the real 
4-dimensional system we are interested in. There is a phase transition from a spin-wave to a 
Kosterlitz-Thouless phase at k ~ 0.5 which is certainly in a different universality class compared 
to the phase transition in four dimensions. In addition there is no spontaneous symmetry 
breaking in the XY-model. However, if on the finite lattice a vacuum expectation value is 
defined through the standard rotation technique |52) . 



«=<^IE^I> ( 105 ) 

then the so defined vacuum expectation value v is zero in the Kosterlitz-Thouless phase but 
v > in the spinwave phase. In this way the 2-dimensional YY-model fakes the situation that 
appears in spontaneous symmetry breaking in 4 dimensions. Increasing the lattice volume, v 
approaches zero also in the spinwave phase as it should. For a finite lattice we expect therefore 
that the fermion mass follows the perturbative relation mp oc yv, if the volume is not too large. 
We will in the following refer to the spinwave phase as the broken and to the Kosterlitz-Thouless 
(or Vortex) phase as the symmetric phase. 

All simulations presented below have been performed in a (2 + l)-dimensional setup with 
a lattice of size L 2 L S with L = 12, L s = 26 being typical values of the lattice extensions. The 
Wilson parameter r has always been set to r = 1 and the domain wall height was chosen to be 
mo = 1, too. Typically several thousand scalar field configurations have been obtained using 
a cluster algorithm. The fermionic quantities have been obtained by standard techniques like 
Conjugate Gradient or Lanczos. In order to avoid convergence problems of these algorithms 
antiperiodic boundary conditions for one of the 2-dimensional indices have been chosen. 

6.1 Mass spectrum 

The fermion spectrum is a first indication for the presence of strong coupling behaviour. 
Fermion masses can be obtained from the propagator 

41 











Figure 10: The fermion mass of the fermion at the waveguide boundary as a function of 
the Yukawa- coupling y. Both, the fermion in the full model (boxes) and the reduced model 
(crosses) are consistent with a linear behaviour, in particular at small y. The dotted line is yv. 
All results are at k = 0.5. 



S st (p) = L- 2 ]T e ip{x ~ y) < >, (106) 



3> y 

xy 



which can be computed using the conjugate gradient algorithm to find matrix elements of the 
inverse fermion matrix on a scalar field background. Using the symmetry 



SUP) = S L L l +2 - s ^ +2 ~\Pp), (107) 

statistics can be enlarged by measuring the LL and the RR components and averaging over 
them. Fermion masses are finally obtained by a fit to a free fermion propagator ansatz 

S{p)rr(ll) = -iZ F (sm(pi) - (+)i sin(p 2 ))/(^sin 2 (p A1 ) +m 2 F ). (108) 



The numerical results can be perfectly described by the ansatz ( |108j ) (see fig. 5 in 

What do we expect for the fermion masses having the phase diagram of fig. 9 in mind? 
Weak coupling behaviour For small values of the Yukawa-coupling mj oc yv with v the vacuum 
expectation value of the scalar field using ( |105| ). The proportionality of the fermion mass to y 



could be confirmed in ||25|| , see fig. 10. In the symmetric phase of the model the fermion mass 
turns out to be zero. This excludes the small y region for a construction of a chiral gauge 
theory as the mirror modes can not be decoupled. 
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Figure 11: The fermion masses from numerical simulations at y = 2 as function of k. The 
boxes correspond to the full and the crosses to the reduced model. 

strong coupling behaviour For larger values of the Yukawa-coupling such that one does not 
enter the region where the large-y expansion of section 5.5 is valid we hope to find our favourite 
PMSi phase. From experience with models showing the strong coupling phase, we expect the 
fermion mass to follow the mean field strong coupling behaviour 

m F oc z~ 2 ; z 2 = Re{V x V* + J . (109) 

The quantity z is non-zero in the symmetric as well as in the broken phase. It is monoton- 
ically decreasing with decreasing k. In the PMSi phase the fermion mass measured from the 
numerical simulation would follow the behaviour in ( |109|) and would therefore increase with 
decreasing k. In fig. 11 the behaviour of the fermion mass as a function of n is shown at y = 2 
where one would expect the strong coupling behaviour to show up if the model would behave 
similar to 2-dimensional Yukawa models with local coupling. 

In the figure, the boxes, connected by a dotted line to guide the eye, correspond to the full 
waveguide model, the crosses to the reduced model. The masses from both models are very 
similar. They stay first almost constant at high values of k and then start to decrease when 
k is lowered. This behaviour is in sharp contrast to the expected increase of the masses if 
there would be a strong coupling region. Unfortunately, the numerical simulations could not 
be driven into the symmetric phase (k < k c — 0.5). The propagators became too noisy and no 
reliable mass value could be obtained. 
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Nevertheless the picture is in clear contradiction with the scenario we have hoped for. The 
mirror fermion masses become small approaching the critical line. This is another indication 
that a PMSi phase does not exist in the waveguide model. Quite in contrary, the fact that 
fermions seem to become light in the symmetric phase at y = 2 point into the direction that 
the fermion spectrum shows weak coupling behaviour. One other important aspect of fig. 10 
and fig. 11 is, that they show the resemblance of the reduced model to the full model. This 
similarity will be relied on when in the next subsection the eigenvalue spectrum is determined, 
which will give an additional hint that the PMSi phase in the waveguide model is missing. 



6.2 Eigenvalue spectrum 

A complementary way to detect a strong coupling phase is the qualitative behaviour of the 
eigenvalue spectrum for different Yukawa models [Q. The shape for a Yukawa-model with 
local (1103) and hypercubical ( |104p couplings are quite characteristic. It would, of course be 
ideal to compute the eigenvalues in the full domain wall model. However, the matrix from 
which the eigenvalues have to be determined is of size 2L S L 2 ® 2L S L 2 . To get all eigenvalues 
of this matrix is almost impossible for reasonably large lattice sizes. We can, however, rely 
on the similarity of the reduced model to the full waveguide model and try to determine the 
eigenvalues in the reduced model at different values of y. 

In fig. 12 the eigenvalue spectra for the reduced domain wall model (a), the model with local 
(b) and with hypercubical (c) Yukawa-coupling is shown. They have been obtained at y — 0.2, 1 
and 4 from left to right and were averaged over 5 scalar field configurations. The lattice size 
has been L — 12. For the waveguide model the eigenvalues at small Yukawa-coupling are very 
similar for all three models and form a narrow band. For y = 1 there appear some differences 
in that for the waveguide and the hypercubical model the band broadens whereas in the model 
with local Yukawa-coupling they seem to develop a gap. 

The difference becomes much more pronounced for y = 4. The local Yukawa-coupling model 
shows a ring structure typical for a strong phase ||63fl . On the other hand it is all too obvious 
that the domain wall model follows the trend of the hypercubical model. Even at y — 4 which 
is a quite strong Yukawa-coupling, no sign of strong coupling behaviour can be seen. The 
eigenvalue distribution just broadens with increasing y. The good representation of the full 
model by the reduced one suggested in figs. 10 and 11 leads to the conclusion that also in the 
full domain wall model the eigenvalues behave very similar. 

This is strengthened by the behaviour of the number of conjugate gradient iterations to reach 
a given residue. This is determined by the condition number of the fermion matrix. In the 
model with a local coupling which has the strong symmetric phase, the number of conjugate 
gradient inversions first rises with increasing y but then, when entering the strong coupling 
region it decreases again, due to the fact that the fermions become massive. No such decrease 
can be seen in the full waveguide model. The number of conjugate gradient inversions just 
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Figure 12: The eigenvalue spectrum for the reduced waveguide model (a), the Yukawa model 
with local (b) and with hypercubical (c) Yukawa-couplings. The values of y are y = 0.2, 1.0 
and 4.0 from the left to the right. 

continues to increase. This is, of course, in full accordance with the results from the large-?/ 
expansion and from the fermion mass behaviour. They do not give a region with strong coupling 
behaviour but in contrast predict light fermions all the way to y — > oo. 

6.3 Dynamical fermions 

All the numerical results discussed in the previous section were obtained in the quenched ap- 
proximation. In the full unquenched model there would have been a single right handed fermion 
coupled to the U(l) gauge field if the fermions would have been in an anomalous representation. 
However, the model in the quenched approximation does not distinguish between an anomalous 
and an anomaly free representation. We could have equally well added an extra fermion species 
S = xM^x to the original action S = ipMip with ip and x having opposite handedness. In the 
quenched approximation the extra fermion x is irrelevant and the quenched approximation of 
the full model with if) and x reduces to the one studied in the previous sections. 

Of course, it would have been desirable to study the waveguide model also dynamically. This 
was also attempted. A Standard Hybrid Monte Carlo algorithm was used in the simulations. 
Using the scalar vacuum expectation value v as an order parameter, the unquenched simulations 
revealed the weak symmetric phase for small values of y, where v = 0. Keeping k fixed, the 
phase diagram in fig. 9 suggests that by increasing y one would enter a symmetry broken phase 
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with v > 0. This pictures is motivated by the ordering effect of the fermions which shift the 
phase transition of the pure bosonic theory towards smaller values of k. Increasing y even 
further, v becomes zero for some critical Yukawa- coupling y c , if one enters a second symmetric 
phase. Whether this is a phase with strong coupling behaviour could then be clearified by 
computing the fermion spectrum. 

This scenario can be tested by numerically computing v using the rotation technique. Un- 
fortunately, for medium size values of y ~ 1, the step size in the Hybrid Monte Carlo algorithm 
had to be decreased by an order of magnitude resulting in large autocorrelation times. This 
became even worse when increasing y further such that no reliable simulation could be per- 
formed in the interesting y region. The bad behaviour of the Hybrid Monte Carlo algorithm is, 
however, consistent with the development of the eigenvalue spectrum discussed in the previous 
section. If a strong symmetric phase exists one would expect the behaviour of the Hybrid Monte 
Carlo algorithm to become much better again since the fermion would have become massive. It 
might be an interesting prospect to attack the difficult region in y again by using the recently 



suggested new simulation algorithm for fermions by Luscher fl65|. 



7 Infinite extra dimension 

As we have seen in the previous sections, our attempts to achieve a chiral gauge theory with 
domain wall fermions seem to fail. The effect of the anti-domain wall can not be switched off 
in case of 5-dimensional gauge fields. In the waveguide model where gauge fields are kept 4- 
dimensional the anti-wall can be shielded from the wall but only for the expense of introducing 
a wave guide wall. New light mirror fermion modes appear on this wall rendering the theory 
vectorlike again. 

We found on the other hand in section 2.1 that on the infinite lattice there is only one chiral 
zeromode and no necessity for introducing a second wall emerges. One might first discard an 
infinite lattice as not leading to a practical proposal for simulations. This is certainly the case 
for the euclidean action formulation of a theory. Narayanan and Neuberger [^, ^j(J suggested 
to take a different route of thinking. They regard the extra dimension as the "time-axis" of a 
(4 + l)-dimensional system and propose to study the 4-dimensional Hamiltonian of this theory. 
The extra dimension is thought of being kept strictly infinite and as such it fulfills the sole 
purpose to project onto the ground state of the system. There will be two Hamiltonians H ± 
corresponding to the two sides of the wall and accordingly there will also be two ground states 
|0± >. The theory can then be defined as the overlap of these two ground states. 

Also in this approach the gauge fields are kept 4-dimensional. As opposed to the waveg- 
uide model there is no need of adding an extra wall and thus no additional fields have to be 
introduced. To keep the extra dimension infinite may also be motivated from the index of the 
fermion matrix M. What one would like to achieve in order to obtain a chiral theory is that 
the fermion matrices M and M\ belonging to the left and righthanded modes, respectively, 
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have different indices. In particular, M should have index 1 and index 0, which would 
correspond to a single chiral fermion as an eigenstate of M. Such a situation can only be 
reached when the matrices M and M< are infinite. Of course, in the domain wall model with 
a finite M it was hoped that one can keep both chiral modes on the lattice such that they can 
be separated and shielded from each other. The model proposed in [25| can be considered as 



a combination of domain wall fermions as suggested by Kaplan [13[] and the idea of infinitely 
many regulator fields advocated by Frolov and Slavnov |)7|]. There it was suggested to intro- 
duce infinitely many Pauli-Villars fermion fields to regulate a chiral gauge theory. Interpreting 
the extra dimension as a flavour space, the different flavours can be associated with new fields 
that have been introduced. Choosing the extra dimension to be infinite, these new fields play 
the role of the infinitely many Pauli-Villars fields in the approach of p7[ . 



To be specific, we will consider the following 4-dimensional action 

X,S /i=l 

+ \ [<W(i + 75) + <w, s (i - 75) - 25s, s] r x • (HO) 



x.s 



This action corresponds to the action (f78|) with the 5-dimensional gauge fields V set to unity 
and the Wilson coupling r set to r = 1. The domain wall mass is chosen as 

m(s) — d — msign(s + -) (HI) 

with d the (odd) dimension of the system. Note that the wall in this setup has been chosen to 
be in the middle of the link connecting the regions in s with positive and negative domain wall 
mass. Of course, keeping the summation over s infinite leads to divergences. However, these 
divergences come from the bulk and one may define an interface effective action 

s t (U) = s eff (u) - i [st ff (u) + s; ff (u)} . (112) 



Sg~ff(U) are the bulk actions for a constant domain wall mass to the right and the left of the 
location of the domain wall. The interface effective action is defined on the wall. It is now 
finite and can be implemented order by order in perturbation theory. 

Staying in the euclidean action formulation would, however, not lead to any practical way for 
numerical simulations. The path Narayanan and Neuberger are following is to construct Si(U) 
from a vacuum overlap formula. This is obtained in a first step by constructing the transfer 
matrix T = e H (with the lattice spacing a = 1) from the lattice Hamiltonian H, H = a)Hd 
with a) and a are the usual fermionic creation and annihilation operators and H is represented 
by a matrix, the explicit form of which will be given below. 

For Wilson fermions with r = 1 the transfer matrix describes the evolution of the system 
from one timeslice to the next and therefore over the smallest distance on the lattice. As the 
transfer matrix contains the full physical information of a system, a description of a model in 
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terms of the transfer matrix is completely equivalent to its path integral formulation. Indeed, 
the path integral Z may be written with the usual trace formula Z = TrT N when N is the 
extension of the lattice in the "time" direction. The construction of the transfer matrix, which 
is a non-trivial task, was performed by Narayanan and Neuberger, following [68[] closely. They 



find the transfer matrices on both sides of the wall, T + and T_ . T + and T_ are the main building 
blocks of the theory in the infinite extra dimension formulation and for the construction of the 
overlap formula. 

As a result a representation of the interface effective action in terms of the transfer matrices 
can be given 

e ^> = iim <i-\(T-rHt+y-^> 



'< b- |(T_) 2 *- 2 |6- >< b+ |(T+) 2 *- 2 |&+ > 

where < 6±| denote some boundary conditions at s — ±00. Choosing these boundary conditions 

as 

|6± >= |0± > (114) 

in the s — > 00 limit the ground state is projected out in ( |113[ ) and the interface effective action 
is given by a simple overlap formula 

e sm =< _ |o+ > . (115) 



To get a feeling for the formula (p. 15 ) which corresponds to the -infinite- lattice effective 



action, we will now as a first step discuss its continuum limit for a (2 + l)-dimensional system. 
We will check in continuum perturbation theory whether we will find the anomaly we expect, 
if the overlap formula correctly describes the chiral gauge theory it was constructed for. The 
Hamiltonians for the left and the right sides of the wall, H- and H + , respectively, are in 
momentum space 

f d 2 k 

- H±(p, q) = H^(p)5 p>q + J [z k T8 qtP+k + z£r t <5, lP+ k| (116) 

where the unperturbed Hamiltonians H are diagonal in momentum space 

fl?(p)=i t" V ?! )- < U7 > 

Pi + ip 2 ±?71 J 

(118) 

\ 1 u 1 

and 

Zk = J d 2 xe ikx [Ai(x) + iA 2 {x)\ . (119) 

If we denote by ^ the set of eigenvectors belonging to positive eigenvalues of the full Hamil- 
tonian 

- H ± *$ = A^J ; Aj>0 (120) 
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we can define an overlap matrix as 



0„ = + . (121) 



P 1 

The continuum interface effective action is then given by 

e s ^ = detO . (122) 

To proceed we will now sketch the perturbative evaluation of Si(A). For this purpose we expand 
the fields 

% = % + i)~C qv + X qD q p (123) 

where ip^ 1 and x are the eigenvectors corresponding to the positive and negative eigenvalues 
of the unperturbed Hamiltonian H^ 1 , respectively. A, B, C and D are coefficients which are 
linear in the gauge fields. The eigenvalues are also expanded 

AJ = A P + AJ« (124) 

with A p the eigenvalues belonging to the free Hamiltonian. The perturbation Hamiltonian can 
be read off from (|116|) 

- Hi(p, q) = z k r8 q>p+k + zlT ] 5 q , p+k (125) 

and the coefficients and the eigenvalues can be obtained from standard Schrodinger perturbation 
theory 

A _ ^Hrtt 
pi - x x Pr? 

A q — A p 



B pq 



\ + A p 



pi ~ ^ \ P^Q 

Aq — Ap 

D P q = \~^ L - (126) 

Ag T Ap 

The diagonal terms A pp and C pp in ( |126| ) are undetermined. This can be partly repaired by 
imposing orthonormality of the eigenvectors which leads to 

A pp + A* p = ; Cpp + C*p = (127) 

which fixes the real part of both coefficients. However, the imaginary part of A and C remains 
undetermined and therefore there is a phase ambiguity in the interface effective action. This 
phase ambiguity reflects the fact that it is the imaginary part of the complex action where the 
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chiral nature is hidden ||69| . The real part of the action is vectorlike. This can also be seen from 
( |113| ). Taking there periodic boundary conditions in the extra s direction, in the limit s — > oo 



< - 10+ > 



'1281 



The effective action ( 128[) is real and corresponds to the domain wall setup with two chiral 
zeromodes of opposite chirality. Thus we see that the imaginary part of the action is the 
crucial ingredient for chiral gauge theories. One proposal to compute the imaginary part of the 
action can be found in |70] . 

The still undetermined phases for the interface effective action have to be fixed. This can be 
done in perturbation theory by demanding that the projection of the true eigenstate onto the 
perturbed one ^^ipp being real. This amounts to Brillouin-Wigner perturbation theory. Of 
course, also for a non-perturbative investigation of a chiral gauge theory the phase ambiguity 
appears. The authors in |2{J propose to take the same phase fixing condition as in perturbation 
theory. The advantage of this condition is that it can be easily implemented and thus the theory 
is completely determined. It is, however, an important -and still open- question whether the 
Brillouin-Wigner condition is suitable also for non-perturbative and strongly fluctuating gauge 
field configurations. 

To arrive eventually at the anomaly one has to go to second order in perturbation theory. 



This can be done |29| and as a result the anomaly is given as 

81 



SA^(x) 



where 



16vr 



E 



Pi + W2 * * Pi 

z z — 

Pi - ip2 



%P2 



-p~-p 



Pi + ip2 



ZpZ—p 



and z is given in ( |119| ) 



A, p + iA 



2,p i 



Ai - - iA 



2.p 



Thus the anomaly becomes 



1 1 

{P2M, P - PiA 2 , p ) = -—F 12 {p) 



(129) 



(130) 



(131) 



(132) 



Eq. (|132|) is the continuum form for the anomaly on the wall from the overlap formula ( 115 ). 
Since in deriving this equation the variation of the action with respect to the gauge fields was 
taken for the effective interface action we should have expected to find the anomaly which 
corresponds to the chiral mode on the wall. This would be the consistent anomaly as we 
discussed in section 4. In fact, the strength of the anomaly in ( |132| ) is 1/47T which is just half of 
the value of the covariant anomaly we have discussed earlier for the domain wall system. Thus 
it corresponds to the consistent anomaly which is the correct result. 
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As a next step it would be certainly interesting to see, whether the overlap formalism can 
be taken over to the lattice. The transfermatrix can be explicitly constructed 

eH±= \C^% C^C C +B±) ■ (133) 
The matrices B and C depend on the gauge fields and their explicit forms are 

B ± = (d + 1 T m) - \ £ \5 X , X+ ^ + 5 X ^U* X _^] (134) 

Z A 4 



c = \ H [8x*>+pU Xt n - 5 X -^U*_^] 7^ . (135) 



2 , 

The effective action on the interface is 

e Si(u) = u < ~ l°+ >u a\<s> + (t/)-*- (u)] / 136 x 
1 <0-|0+> 1 ' V ; 

The |0± >u are the ground states of the Hamiltonian in presence of the gauge fields whereas 
|0± >i denote the ground states for free fermions. The phases Q±(U) are undetermined which 
again reflects the phase ambiguity we have encountered in ( [126|J127| ). The phases are fixed 
using the same condition as in perturbation theory also for the non-perturbative case 

c » ±(t = ^ <0± I Q± >i . (137) 
| [7 <0±|0±>i| v ' 

In close analogy to the earlier discussion of the anomaly an external gauge field is chosen 

U(k)=e^ c °<^ n+ ^ (138) 

with 

A^) =A M + 20sin(^) . (139) 

The anomaly equation now reads 

_« <^ i 

7T L^oo dq> A-k 

In order to compute the anomaly numerically, the Dirac sea has to be filled and the eigenfunc- 
tions corresponding to the negative energy eigenstates have to be found. This amounts to find 
the eigenvalues of ( |133| ) numerically. In the computation was taken to be A2 = 0.32. In P^j 
two values of the domain wall mass were chosen, m = 0.5 and m = 0.9 and the anomaly was 
computed on different size lattices. At both m the value of the anomaly was extrapolated to 
L = 00. Both extrapolated values coincide and give —0.02545(5). This is to be compared to 
the value of the anomaly in the continuum of —0.02546. The same test of the anomaly equation 
to be fulfilled on the finite lattice can also be performed in four dimensions by a clever choice 
of the external gauge fields keeping the Hamiltonian to be at least block diagonal. Again the 

51 



continuum value of the strength of the anomaly could be reproduced in the thermodynamic 
limit. 

We discussed in this section a proposal of keeping the extra dimension in the domain wall 
model strictly infinite. Alternatively we might think of having introduced infinitely many regu- 
lator fields for a chiral gauge theory. We saw, that in continuum perturbation theory the correct 
value of the consistent anomaly is found. Even more remarkable is that the overlap formula 
resulting from the ground state projection can be implemented on the finite lattice. The task 
there is to diagonalize a -admittedly huge- matrix. Choosing external gauge fields that vary 
smoothly, the value of the consistent anomaly could be reproduced when the finite lattice re- 
sults are extrapolated to the infinite lattice. With suitable techniques of matrix diagonalization 
a numerical simulation may be reached eventually. 

With this outlook it might be possible to test whether a chiral gauge theory can be obtained 
when the gauge fields are made dynamical. Of course, the tests being performed on the overlap 
formula so far correspond to the situation in the domain wall model where everything worked. 
There we also found for smooth external gauge fields that the correct -in this case covariant 
anomaly- can be reproduced on the lattice and that the chiral zeromode spectrum comes out 
as desired. This encouraging picture was only destroyed when dynamics had been added. 
Therefore it seems that only a full simulation can decide the question whether the infinite s 
approach is the solution of the chiral fermion problem on the lattice. 

In a numerical simulation only the real part would be used for the updating. Given a new 
gauge field configuration consisting of a change of only one link, the imaginary part would be 
computed by diagonalizing H±. It is then hoped that, since only one link is changed at a time, 
the old eigenfunction of H± is a good starting point for the new diagonalization process. The 
general idea is not to take the imaginary part in the updating process but into the measured 
observables. It is, of course, not clear at all that the fluctuations of the imaginary part for 
realistic lattice gauge field configurations are small. If these fluctuations are large, which might 
well be possible due to the roughness of the 4-dimensional gauge fields, this would not lead 
to any practical simulation algorithm. On the other hand an implementation of the above 
simulation proposal could provide the first measure for the strength of the fluctuations of the 
imaginary part of chiral fermion actions. 

One might find the idea unsatisfactory that a chiral theory can only be constructed by 
considering the Hamiltonians for both sides of the wall and construct the overlap of the corre- 
sponding ground states. There seems to be at the moment no way of connecting the overlap 
setup with any euclidean lattice formulation of the same theory. However, the anomaly in 
the infinite dimensional lattice approach for chiral lattice fermions can be understood from 



a limiting procedure of a lattice version with finite extend in the s-direction [45|. In this 



version the s-direction is split into a gauged region for s < N and an ungauged region for 
iV < s < N < oo. The limiting procedure which gives the correct -consistent- anomaly is to 
let both N and N' go to infinity simultaneously. Therefore the interface effective action ( |112|) 
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might be understood as the limit s —>■ oo of a gauge variant finite lattice version. As we have 
seen in the discussion of the waveguide model, section 5.2, a gauge variant action can be made 
gauge invariant by adding scalar or Stuckelberg fields. But then the infinite s effective action 
is hiding that it came from the waveguide model on the finite lattice. As the waveguide model 
failed in constructing a chiral gauge theory, the same faith may await the overlap approach to 
chiral fermions on the lattice. 



8 Domain wall fermions and QCD 

Even if the domain wall fermion approach to chiral gauge theories on the lattice may fail, there 
could be a quite unexpected application of this idea. As put forward by Shamir |43|, [H]] domain 
wall fermions could be used for simulating QCD. This seems to be a natural approach since 
the odd dimensional theory one starts with is vectorlike. Consider a slightly modified action 
for domain wall fermions where we choose 2N points in the s-direction. 

s F = EE^P)-a(!/) + H^ 

X,S fl=l 

+ \ E^x [<W 5 (1 + 7s) + - 75) - 25s, s] r x 

z x,s 

+ r#f (1 + 75)^ + ^(1-75)^ • (141) 

This action is identical to eq. (|110|) except that the boundary term in the s-direction is made 
explicit. The gauge fields are kept 4-dimensional. The boundary terms are equipped with an 
additional new parameter fn. For fn = 1, which corresponds to periodic boundary conditions, 
no zeromodes appear. For m = we will have open boundary conditions and the action (|141|) 
describes boundary fermions with chiral surface modes pSfl . 

For a non-zero value of fn an interaction of the surface modes at s = 1 and s = 2N is 
induced and they will combine to a Dirac state. In order to derive the mass dependence of this 
state we compute the propagator from the action ( |141| ). First there is a homogeneous solution 
of the Dirac equation 

G (s,s') =Be- als - s ' 1 (142) 



with 



B' 1 = 26sinha 

1 + b 2 (p) +p 2 



2 cosh a (p) - 

b(p) 

p = sinp fl b(p) = 1 - m + ^(1 - cosQo M )) . (143) 

The full propagator is given as 

G(s, s') = G (s, s') + +„e~ Q(s+s,) + A + e- ai2N ~ s ~^ 

+ A m (e-^ N+s ~ sl) + e -( N +°'~^ (144) 
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where 



A 



A 



A 



in 



+ 



A- 1 J B6(l-m 2 )(e a -6) 
/\- l Bb{l-m 2 )(e- a -I 
2A _1 i?6 2 mcosh a 



A 



be a - 1 + m 2 b(e a - b) . 



(145) 



This complicated expression can be evaluated in the limit of small momenta and small 
values of rft. As a result one finds that the propagator in momentum space is 



For the boundary fermion model m = 0. We see that in the low momentum limit the propagator 
describes a massless particle. The s-dependence of the full propagator ( |144j ) shows that the 
chiral zeromodes appear as surface modes on the boundary of the lattice in the s-direction. 
This is the same result as we discussed in section 2 where the zeromodes were computed in the 
Hamiltonian language. They are depicted in figure lb. 

Turning on TFT, these zeromodes are coupled through the links connecting the last and the 
first s-slice and acquire a -Dirac- mass. As can be read off from ( |146| ) this mass is proportional 
to rft which might therefore be interpreted as an external quark mass. The reason, why this 
result makes boundary fermions interesting for QCD simulations is the following. According 
to the nogo theorems chiral symmetry is broken on the lattice. However, as weak coupling 
perturbation theory shows, chiral symmetry will be restored in the continuum limit of lattice 
QCD. To reach this limit a careful tuning of the bare parameters in the lattice QCD action is 
necessary to approach the phase transition. The tuning of the pion mass to zero as a function of 
the bare parameters is non-trivial and demanding. In particular, there are 0(l/a) corrections 
to the renormalized quark mass. 

Eq.( |146| ) offers a nice alternative. It suggests that quark masses only get multiplicative 
renormalization and that one has therefore only to tune fft to zero in order to reach the phase 
transition. In this context, fft might be thought of as an external source. In particular, the 
phase transition point is expected to be fn = 0. For a finite lattice in the s-direction there 
is an additional complication. Since the wavefunctions of both modes on the walls fall off 
exponentially, there will be always an overlap of both quark modes in the middle of the lattice. 
The value of the quark's wavefunction can be estimated to be (1 — m) s . Therefore the extension 
of the lattice in the extra dimension should be N « —log(m). Whether the advantage to make 
tuning trivial will merit the introduction of an extra dimension can only be answered in actual 
simulations. 

For a finite system the overlap of the surface modes lead to some anomalous effects in the 
middle of the lattice. The currents are defined as 




(146) 




(147) 
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for pi — 1, 4 and 



J5 ~< C^-^r .=2* • (148) 

Since the 5-dimensional theory is vectorlike, it is anomaly free and the divergence equation is 
®n3n = ~9 5 j 5 , n = 1, ...,4. Written out explicitly 



— mjl~ 2N s = 1 



■l 



E^= "Otf! Ks<2iV . (149) 

M=l | _ 3 s=2N-l _ mf =2N g = 2N 



One can define an axial current as in section 3 

2N 



A M = -]Tsign(iV- S + ^. (150) 

s=l ^ 



Then the divergence is 



= 2mj! N + 2jf . (151) 
Whereas appears to be derived from a usual classical massterm, j% N is an additional anoma- 



lous term stemming from the overlap of the surface modes. As is shown in |71j this term vanishes 
with iV — > oo as 

exp(-c /(7 2 )exp(-AiV) (152) 

with Co a constant, go the gauge coupling and A the confinement length. From this formula 
it is seen that for every g Q fixed, the anomalous term vanishes if the extent of the lattice in 
the extra dimension is sent to infinity. This suggests that one might try to follow the program 
of the previous section and work in the N = oo limit directly. Since in the application of the 
domain wall model to QCD the effective action on the boundary is real and therefore does 
not possess the phase ambiguity we found for the chiral gauge theory case, the program of 
Narayanan and Neuberger can now be performed without subtleties involved. The derivation 
of the overlap formula for the QCD case can be worked out following closely and an explicit 



formula is found in [7T|. The vanishing of the anomalous term in the N — >• oo limit constitutes 



a proof that axial symmetry is restored at the QCD phase transition. In the derivation of this 
result basically no approximation is involved. This amounts to the statement that restoration 
of axial symmetry will also appear at strong coupling and is valid not only to every order in 
weak coupling perturbation theory |l|] but holds also non-perturbatively. 

9 Summary and conclusions 

In this article we summarized the status of domain wall fermions invented to construct a chiral 
gauge theory on the lattice. In the continuum domain wall fermions are obtained as chiral 
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zeromodes bound to a -soliton shaped- mass defect living in an extra odd dimension. The 
domain wall model can be taken over to the finite lattice by choosing the mass defect to be a 
step function. Through the necessary introduction of boundary conditions a second anti-wall 
is generated and the chiral zeromodes appear pairwise as left and righthanded fermions. They 
are bound to the wall and the anti-wall, respectively, and fall off exponentially with increasing 
distance from the walls. Alternatively one may choose free boundary conditions in the extra 
dimension in which case the chiral modes appear as surface modes. 

The charge flow in the domain wall model is according to the picture developed in the 
continuum by Callan and Harvey. The anomaly appearing in the zeromode current along 
the wall is compensated by the divergence of the Chern-Simons current induced by the heavy 
fermions that live off the wall in the extra dimension. The Callan-Harvey analysis can be 
performed on the lattice confirming the continuum picture at least, as long as the gauge fields 
are taken to be smooth and external. The correct strength of the -covariant- anomaly can be 
directly measured on the finite lattice. 

The crucial step for the realization of a chiral gauge theory on the lattice is the coupling 
to dynamical gauge fields. The fluctuations of the gauge field can become strong and the fear 
is that they will lead to an interaction of the two chiral zeromodes combining them to a Dirac 
particle. In the waveguide model only one of the two domain walls is gauged. Interpreting the 
extra dimension as a sophisticated flavour space, the gauge fields are kept strictly 4-dimensional 
and the same for all flavours in a region around the "wall" -flavour, i.e. the location of the wall. 
At the boundaries of this "waveguide region" gauge invariance is lost. This can be repaired by 
introducing scalar or Stiickelberg fields at the boundary of the waveguide region. The coupling 
of the scalar fields and the fermions are of the Yukawa-form and will be equipped with a Yukawa- 
coupling y. For y = one essentially obtains the model with free boundary conditions and 
new fermionic zeromodes appear which can interact with the domain wall zeromode through 
the gauge fields. The important and crucial question for the model to succeed is then whether 
the fermionic spectrum at the waveguide boundary can be made massive such that the modes 
decouple in the continuum limit and one is left with only the chiral zeromode on the gauged 
wall. 

The result of the investigations in the domain wall model can be given, comparing fig. 9 and 
fig. 13. These figures show the phase diagrams that would have led to a possible construction 
of a chiral gauge theory with domain wall fermions (fig. 9) and the much more realistic phase 
diagram in fig. 13 which prohibit such a construction. The phase diagram in fig. 13 was obtained 
by a large-?/ expansion in combination with numerical simulation results. 

The interesting part of the phase diagram in fig. 9 is the PMSi phase. In this phase the 
fermions at the boundary of the waveguide region would have become massive. There the 

3 In Fig. 13 the PMS2 phase is left out. The existence of this phase was not investigated. But even if it would 
exist, a construction of a chiral gauge theory would not be possible because all modes including the one on the 
domain wall become heavy. 
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Figure 13: A realistic phase diagram for domain wall fermions. There will be a symmetry bro- 
ken (FM) phase as well as two symmetric phases (PMWi and PMW2) with different vectorlike 
spectra of massless fermions. 



Yukawa-coupling would have been strong enough to bind the fermion and the scalar field to 
form a massive bound state. Therefore this state would have decoupled from the low energy 
spectrum, leaving behind a chiral zeromode on a domain wall coupled to gauge fields in a gauge 
invariant manner. In earlier investigations of coupled fermion-Higgs models such a phase could 
be established or it was missing, depending on the particular implementation of the form of the 
Yukawa-coupling. What is found in the waveguide model from the numerical simulations and 
the large-?/ expansion is that the PMSi phase is absent. Instead, there are two symmetric phases 
with weak coupling behaviour both of them with a vectorlike spectrum of zeromodes. Thus in 
neither of these symmetric phases a chiral gauge theory can be obtained. A word of care has 
to be given at this point. Although the phase diagram in fig. 13 represents a consistent picture 
from the large y expansion and numerical simulations, the simulations were very difficult in 
the most interesting region of the phase diagram. Although evidences are strong that a PMSi 
phase is absent, its non-existence could not be proven without any doubt. Maybe using the 
algorithm proposed in |)5 can help to achieve better simulation results in the difficult region. 

What remains from the domain wall model idea is a variant where the extra dimension 
is kept strictly infinite. In this interesting line of thinking the extra dimension is regarded 
as the time axis of a 4-dimensional Hamiltonian. The Hamiltonian is split into two parts 
corresponding to the two sides of the domain wall. Due to the infinite extension in the extra 
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dimension the ground states on both sides of the wall are projected out. The resulting theory 
consists of the overlap of these two ground states. Assuming smooth external gauge fields it 
could be shown in perturbation theory and non-perturbatively through a diagonalization of 
the finite lattice Hamiltonians that this formulation of the domain wall model gives again the 
correct -consistent- anomalies. So far, in this interesting approach no dynamical gauge fields 
have been considered and it remains to be seen whether the promising results survive their 
inclusion. 

A completely different application results from the interpretation of the wall zeromodes as 
massless quarks. The boundaries of the lattice in the extra dimension can be coupled through 
a mass term. In this way, the surface modes combine to build a Dirac fermion. The interesting 
aspect of this approach is that the mass of the Dirac-particle is directly proportional to the 
new mass parameter. Thinking of simulations in QCD this would mean that one could perhaps 
avoid the fine tuning problem of the quark mass in the conventional formulations of lattice 
QCD. 

To conclude, the clash between chiral gauge invariance and the existence of a single Weyl 
fermion on the lattice could also not be resolved by the domain wall fermion approach. It seems 
that non-perturbative lattice models of chiral gauge theories unavoidingly lead to doubler states 
which might be interpreted as mirror fermions. One is therefore left with the choice to believe 
in some -yet undiscovered- new fermions, or to think further, maybe along the lines of the 
Rome-approach JT^], how nature managed to construct chiral gauge theories. 
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